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Abstract. Entanglement in J\ — J2, S = 1/2 quantum spin chains with an impurity 
is studied using analytic methods as well as large scale numerical density matrix 
renormalization group methods. The entanglement is investigated in terms of the von 
Neumann entropy, S = — Ty pa logpA, for a sub-system A of size r of the chain. The 
impurity contribution to the uniform part of the entanglement entropy, Si mp , is defined 
and analyzed in detail in both the gapless, J 2 < J$i as wen as the dimerized phase, 
J2 > Jf , of the model. This quantum impurity model is in the universality class of the 
single channel Kondo model and it is shown that in a quite universal way the presence 
of the impurity in the gapless phase, J 2 < J2, gives rise to a large length scale, 
associated with the screening of the impurity, the size of the Kondo screening cloud. 
The universality of Kondo physics then implies scaling of the form Si mp {r/^K, r / R) f° r 
a system of size R. Numerical results are presented clearly demonstrating this scaling. 
At the critical point, J%, an analytic approach based on a Fermi liquid picture, valid at 
distances r ^> £k and energy scales T -C Tk, is developed and analytic results at T = 
are obtained showing Si mp = + n(l — r/R) cot(7rr/i?)]/(12i?) for finite R. For 

T > 0, in the thermodynamic limit, we find Si„ lp = [ir 2 £kT / (6v)} coth(27rrT/w). In 
the dimerized phase an appealing picture of the entanglement is developed in terms of 
a thin soliton (TS) ansatz and the notions of impurity valence bonds (IVB) and single 
particle entanglement (SPE) are introduced. The TS-ansatz permits a variational 
calculation of the complete entanglement in the dimerized phase that appears to be 
exact in the thermodynamic limit at the Majumdar- Ghosh point, J2 = Ji/2, and 
surprisingly precise even close to the critical point J%. In appendices the TS-ansatz is 
further used to calculate (S*) and (S r ■ SV+i) with high precision at the Majumdar- 
Ghosh point and the relation between the finite temperature entanglement entropy, 
S(T), and the thermal entropy, Sth{T), is discussed. Finally, the alternating part of 
Si m p is discussed, together with its relation to the boundary induced dimerization. 
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1. Introduction 

Much of the mystery and power of quantum mechanics arises from entanglement, which 
leads to Einstein's "spooky action at a distance" but is now recognized as a resource 
by the quantum information community, being essential for quantum teleportation 
or quantum computing [1]. Ground states of quantum field theories and many 
body theories exhibit fascinating entanglement properties which are beginning to be 
understood [21 El IU El El E] • A useful measure of many body entanglement when the 
total system is in a pure state is the von Neumann entanglement entropy This is 
obtained by focusing on bipartite system where space can be divided into 2 regions, A 
and B. Beginning with the ground state pure density matrix, region B is traced over 
to define the reduced density matrix p A . From this the von Neumann entanglement 
entropy [HI E] , 

S{r) = -Tr\p A hip A ] (1.1) 

is obtained for subsystem of size r. Several other measures of entanglement are in 
current use such as the concurrence [TOl [TT] which is monotonically related to the 
entanglement of formation [12], the distillable entanglement [121 03] and the relative 
entropy of entanglement [TH [151 HE]- See also p2]. Some of these measures have been 
developed to describe entanglement as it occurs in systems in a mixed state where it is 
much harder to quantify entanglement, see for instance the seminal paper by Bennett 
et al. [12]. Here we focus mainly on bipartite systems in pure states for which the von 
Neumann entropy, S, is an essentially unique measure of the entanglement. The rate at 
which S grows with the spatial extent, r, of region A is not only a fundamental measure 
of entanglement, it is also crucial [HI [191 EQl 121] to the functioning of the Density 
Matrix Renormalization Group (DMRG) a powerful numerical method for solving many 
body problems [221 123]. For systems with finite correlation lengths, it is generally 
expected that S A grows with the area of the boundary of region A [2U ES]- In the one- 
dimensional case, conformally invariant systems (with infinite correlation length) have 
S(r) — > (c/3)lnr [21(26] where c is the "central charge" characterizing the conformal 
field theory (CFT). Entanglement entropy has recently been shown to be a useful way 
of characterizing topological phases of many body theories [23 [23 (2H1 [30] , which cannot 
be characterized by any standard order parameter. Entanglement entropy is also closely 
related to the thermodynamic entropy of black holes and to the "holographic principle" 
relating bulk to boundary field and string theories [23 EH E2] • Hence, due to these latter 
developments, even though the von Neumann entanglement entropy may only give an 
incomplete description of entanglement in mixed states, such as would be the case at 
finite temperature or in the presence of noise from the environment, an understanding 
of its behavior in such states is important from the condensed matter perspective. We 
discuss the precise connection between the finite temperature entanglement entropy and 



the thermal entropy in some detail in Appendix E 



Recent experiments [33] on the magnetic salt LiHo x Y 1 _ a; F4 have been interpreted 
as evidence for quantum entanglement in the magnetic susceptibility and electronic 
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specific heat at temperatures approaching IK. Theoretical work have shown that 
macroscopic entanglement is in principle observable at much higher temperatures [31] 
and should also be observable using other probes such as neutron scattering [35J. 
Experimental measures of the concurrence have thus been obtained [33l [35] . Some 
of the theoretical [361 EZl EHl EHl SOI [HJ S21 S3] and experimental [H 05] have focused 
on establishing entanglement witnesses (EW) for detecting entanglement. An EW is 
an hermitian operator W with (p)w — tv(Wp) > for all separable p. Thus, if 
(p)w < 0, then p is non-separable. Energy fluctuations [39], the persistent current [39|, 
the temperature [13] and the magnetic susceptibility [HH 02j [35] have been proposed as 
EW's. While an EW can detect entanglement it does not provide a quantitative measure 
of the strength of the entanglement. Analysis based on the magnetic susceptibility as 
an EW have been interpreted as a signature of entanglement at temperatures as high 
as 365K in the nanotubular system Na2Vas07 [45], close to 100K in the warwickite 
MgTiOB0 3 [B] and around 20K in the pyroborate MgMnB 2 6 [B]. 

Comparatively few results [26j H7J HH1 SHI EH EH E2, [53] have been obtained on 
entanglement in systems with impurities. For CFT's impurity interactions will, under 
the renormalization group (RG), flow to fixed points which can be represented by 
conformally invariant boundary conditions. These conformal boundary conditions can 
be characterized by the zero temperature thermodynamic impurity entropy [51], hig, 
a length independent term in the thermodynamic entropy of a semi-infinite system. 
Consequently, it was argued [26] that the entanglement entropy of a semi-infinite CFT 
is (c/6) ln(r/a) + lng, where a is a constant which is non-universal but independent of 
the boundary condition. At the fixed point, the \ng term is then part of the impurity 
contribution to the entanglement entropy. The RG flow between different fixed points 
was numerically studied in terms of In g for the quantum Ising and XX Z chain in |47j . 
It is then of considerable interest to quantitatively define what we call the "impurity 
entanglement entropy" : 

Simp = S^with impurity) — S^no impurity), (1.2) 

the additional entanglement entropy that arises from adding an impurity in region A. 
A specific implementation of Eq. (11.21) suitable for our numerical work will be discussed 
in section [21 

Here we focus on the impurity entanglement entropy in the Kondo model and a 
related spin-chain model. The 3-dimensional (3D) Kondo model Hamiltonian is: 

H= [ d 3 r[^{-V 2 /2m)i> + J K 5 3 (r)^(a/2)ij ■ S], (1.3) 



where tp(r) is the electron annihilation operator (with spin-index suppressed) and the 
S a are S — 1/2 spin operators. Most of our numerical work is performed on the 
related J\ — J2 family of antiferromagnetic S = 1/2 Heisenberg spin chain models. 
For J 2 < J| — 0.2412 the J\ — J2 spin chain is gapless and the low energy behavior of 



an impurity is equivalent to that of the Kondo model. In Appendix A we review the 



connection between these models. (This connection is pursued further in [55].) In the 
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J\ — J 2 spin chain model the equivalence to the Kondo model is achieved by modeling 
the impurity as a weakened coupling, J' K , at the end of an open chain, (see Fig. [TJ) We 
then write the hamiltonian for the spin chains as: 

R-l R-2 

H = J' K [S x ■ S 2 + J 2 S X ■ 5 3 ) + Y,Sr ■ S r+1 + J 2 ^SV ■ S r+2 . (1.4) 

r=2 r=2 




Figure 1. Schematic picture for the J\ — J2 spin chain model (|1.4p with an impurity 
located at the left boundary and coupled with J' K . 

While unimportant for our analytic work, the spin chain representation, Eq. (jl.4p . 
dramatically alleviates the numerical work compared to using the Hamiltonian, 



Eq. (11.31) . As discussed in Appendix A, the finite size corrections are much smaller 
when J 2 is fine-tuned to the critical point J c ~ .2412 [56], a fact that, from a numerical 
perspective, presents a considerable advantage. For J 2 > J 2 the spin chain enters a 
dimerized phase [57] with a gap and the relation between Eq. (11.41) and Kondo physics 
no longer holds. The two-fold degenerate dimerized ground-state is exactly known at the 
Majumdar- Ghosh [58] (MG) point, J 2 = J/2, with the two ground-states corresponding 
to the two possible dimerization patterns. In the dimerized phase the fundamental 
excitations can be viewed as single 5 = 1/2 solitons separating regions with these two 
distinct dimerization patterns. A sketch of the phase-diagram of the J\ — J 2 model 
summarizing these points is shown in Fig. [2j 

At low energies, the spin chain model is equivalent to the Kondo model in any 



dimension, D. This is also true for the entanglement entropy as we show in | Appendix B 
A physical motivation for studying the Kondo model in the D = 2-dimensional case is 
provided by the possibility of using the spins of gated semi-conductor quantum dots 
as qubits [591 EH EI]. The various quantum dot spins would ideally be entangled 
only with each other. However, in practice, they would also be entangled with the 
surrounding conduction electrons which would cause limitations on the functioning of 
such a quantum computer [62], [53]. More generally, any physical realization of a quantum 
computer has a dissipative environment which is entangled, to some extent with the 
qubits. The Caldeira-Leggett (spin-boson) model [63] of a spin interacting with an 
ohmic dissipative environment is equivalent, at low energies to the Kondo model and 
entanglement in the Caldeira-Leggett model has been considered in recent work [62| 164"] . 
The understanding of entanglement as it occurs in these quantum impurity models is 
therefore of considerable interest. 
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Figure 2. Phase diagram of the J\ — Jn model as a function of J2. The critical point 
J 2 ~ .2412 separates a gaplcss Hciscnbcrg phase from the dimerized phase. At the 
Majumdar-Ghosh (MG) point, J 2 = Ji/2, the two exact ground-states, A and B, for a 
system with an even number of spins are shown. For an odd number of spins a single 
spin (soliton) will separate regions with these two dimerization patterns. 



Spin chains have also been proposed for the purpose of creating quantum 
communication channels [65], [66], [67] and subsequent work has showed that entanglement 
in S = 1/2 spin chain models can be used to establish perfect state transfer [681 1551 170] . 
In most cases the couplings in the bulk of the chain are not uniform but vary with r, 
however, a model closely related to Eq. (11.41) . where only the couplings at the end of the 
chain are modified have also been shown to accommodate perfect state transfer |71j . 
The algorithm for perfect state transfer proposed in [68J has been experimentally 
realized using a three-qubit liquid NMR quantum computer simulating a Heisenberg XY 
interaction [72J. The possibility for using quantum spin chains as perfect quantum state 
mirrors has also recently been emphasized [731 [71]. From this perspective, entanglement 
in the spin chain model, Eq (11.41) . is clearly of interest even aside from the relation to 
Kondo physics. 

In the ground state of the Kondo model, the spin of the impurity is "screened" 
meaning that it effectively forms a singlet with the conduction electrons. This screening 
is expected to take place at an exponentially large length scale 



Here v is the velocity of the low energy excitations (fermions or spin-waves) and 
Tk, the Kondo temperature, is the energy scale at which the renormalized Kondo 
coupling becomes large. One of the interesting features of entanglement entropy in 
many body systems is its universality, which makes it useful for studying quantum 
phase transitions [3j HI [5j [26] . A major conclusion of our work is that this universality 
extends to Si mp in the following sense. When is large compared to all microscopic 
length scales in the system, and the system size R is 00, we find that S imp (r/^ K ) is a 



£ K = v/T K oc e 
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universal scaling function depending only on the ratio t/^k- For finite systems of extent 
R we find when r, R — r, ^> 1 : 

Simp = S im p(r/R,r/^ K ) (1.6) 

In this case, there are actually two different scaling functions depending on whether 
the total spin of the finite system ground state is or 1/2. (For the spin chain these 
two cases correspond to R even or odd respectively.) We emphasize that this holds 
for entanglement in both models Eq. (jl.4p and (jl.3p . This extends to entanglement 
entropy the well-known universality of other properties of the Kondo model. We expect 
this universal scaling property to hold generally for quantum impurity models as was 
remarked on in [47]. An immediate consequence of this universality is that the same 
impurity entanglement entropy occurs, for large r, in many different microscopic models. 
It is insensitive, for example, to the dimensionality of space and the range of the Kondo 
interaction (as long as it is finite). In [75] preliminary density matrix group (DMRG) 
results obtained for the spin chain model, Eq. (jl.4p . at J2 = J 2 showed that the r- 
dependence of Si mp confirms this picture and the presence of the length scale was 
demonstrated. In section [2] we provide additional evidence supporting this scaling at 
J 2 = J 2 as we h as results for for J 2 < Jf . 

The universal scaling functions are generally not amenable to analytic calculation 
except in certain limiting cases. The most straightforward of these is when <C r, 
in which case it is possible to perform calculations using Nozieres local Fermi liquid 
theory (FLT) approach [76] . as developed in [771 [78], [79] . in EE] initial analytical results 
using this approach were reported and in section [5] we present a detailed derivation and 
additional results notably at finite temperature. 

When the Kondo coupling approaches either the weak coupling fixed point, J' K = 0, 
or the strong coupling fixed point, J' K — 1, one might have assumed that the impurity 
entanglement entropy would vanish. This turns out to be the case at the strong coupling 
fixed point, J' K = 1. However, as we initially reported in [75], the impurity entanglement 
at the weak coupling fixed point is non-zero. In sections 121171 we discuss in detail this fixed 
point entanglement entropy and further develop the intuitive picture for understanding 
it. 

For J 2 > J| the spin chain enters a dimerized phase [57] with a gap and the relation 
between Eq. (11.41) and Kondo physics no longer holds. However, entanglement as it 
occurs in spin chains without impurities, viewed as model systems for entanglement, 
is currently the focus of intense studies and is relatively well established both from a 
static [HQlIHIlElIHaiHalESlIHSESlIHSlIHIland dynamic perspective [HHl ESI ED] . Here 
we show that it is possible to obtain almost exact analytical results for the fixed point 
entanglement for the J\ — J2 model, Eq. (jl.4p . in the dimerized phase. This approach 
is based on describing the lowest lying excitations in the dimerized phase as simple 
domain walls or single site solitons [91] which we refer to as thin solitons (TS) and is 
described in detail in section [61 While the use of gapped spin chains for the purpose 
of quantum communication and quantum computing is less evident than for gapless 
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chains, the relative simplicity of the entanglement as it occurs in the dimerized phase 
allows for the development of an appealing intuitive picture of how the entanglement 
arises in terms of single particle entanglement (SPE) and impurity valence bonds (IVB). 
These concepts can be rigorously established using the TS approach in the dimerized 
phase and, more importantly, appear to be quite general concepts applicable to other 
models even in the absence of a gap. (See section [2]). 

The outline of the paper is as follows: In Section [2] we discuss in detail our definition 
of Si m p and show additional evidence for the scaling behavior of Eq. ( 11. 6ft at J2 = J%- 
The intuitive picture for understanding the impurity entanglement in terms of SPE and 
IVB is developed in Section [3] along with the fixed point entanglement. Weak scaling 
violations, related to another single-site measure for the impurity entanglement, Sj mp , 
are discussed in Section HI Section [5] contains a detailed account of the FLT approach 
for calculating Si mp and Section [6] describe the thin soliton approach to performing 
variational calculations for the entanglement in the dimerized phase for J 2 > Jf. 
Numerical results for the fixed point entanglement are presented in section [71 Most 
of the numerical results presented are density matrix renormalization group (DMRG) 
calculations performed on parallel SHARCnet computers keeping m = 256 states. As 
we discuss in section [2] our working definition of Si mp only focus on the uniform part of 
the entanglement. In section [8] we therefore present results for the alternating part of 
the entanglement as well as the dimerization. Finally we briefly summarize our main 
results in section [91 

This paper contains a number of appendices, some of which may be of quite general 
interest. In Appendix A we review field theory results on the .D-dimensional Kondo 
model and the spin chain, and their relationship to each other. In Appendix B we extend 
these results to prove that the impurity entanglement entropy is the same for the D- 
dimensional Kondo model and the spin chain model. This appendix also contains a new 



derivation of the free fermion entanglement entropy in D-dimensions. In Appendix C 



we derive the 7-point formula used in the numerical work for extracting the uniform 
and alternating part of the entanglement entropy. In Appendix D we prove that the 
entanglement entropy is the same for all linear combinations of spin up and spin down 
elements of a doublet state. The connections between finite temperature entanglement 
entropy and thermal entropy are discussed in Appendix E In Appendix F we present 
new results on the Majumdar-Ghosh model based on the thin soliton ansatz, including 
(S*) in the ground state with open boundary conditions and an odd number of sites, 
and the dimerization, (S r ■ S r+ \). 



2. The impurity Entanglement Entropy 

We begin by discussing our definition of the impurity entanglement entropy, Eq. (jl.2p . 
For impurity problems such as the Kondo model it is quite standard in experimental 
situations to define the impurity contribution to, for example, the susceptibility or 
the specific heat by subtracting reference values with the impurity absent. We have 
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therefore defined the impurity entanglement analogously. For the S = 1/2 spin chain 
the entanglement entropy has not only a uniform but also a staggered part [87]. For 
r> 1 we can write: 

S(r,R) = Su(J' K ,r,R) + (-l) r S A (J' K , r, R), (2.1) 

where Su{J'ki r i R) an d Sa{J'ki r i R) are slowly varying functions or r. The entanglement 
entropy is also strongly dependent on whether the total spin of the ground-state of the 
system is or 1/2, or, equivalently, whether R is even or odd. 

[In general, for R odd, the ground state is a spin doublet. The entanglement entropy 
is the same for the spin up or down element of the doublet or for any linear combination 
of these states. We give a formal proof of this in Appendix D The result follows from 
the fact that any linear combination of spin up and down can be obtained by a spin 
rotation from the spin up state. It is intuitively obvious that the entanglement entropy 
doesn't depend on the direction of the spin quantization axis. For higher spin states the 
situation is more complex and in general different elements of the spin multiplet have 
different entanglement entropy] 

Theoretical work [21 [26] has established that for a uniform S = 1/2 spin chain with 
periodic boundary conditions, or in general for gapless 1-D models, the entanglement 
entropy is given by: 



S PBC (r,R) = ~\n 
3 



R . fnr\ 
— sin — 
71 \RJ 



si, (2.2) 



with c = 1 and Si ~ 0.726 [92]. In this case there is no alternating term in the 
entanglement entropy. For the purpose of studying quantum impurity models, our 
focus is here exclusively on the case of an open S = 1/2 spin chain where the alternating 
term, Sa, is non-zero [ST] and it is known [26] that for a gapless system of linear extent 
R, the uniform part has a bulk part, 



S uo (r, R) « - In 
6 



~2R 


Cixr \ 


— sm 


\R ) 


71 





\ng + Sl /2, (2.3) 



independent of J' K , as well as an impurity contribution that will depend on J' K and 
that is our main focus. Here lng is the zero temperature thermodynamic impurity 
entropy [51]. We see that the presence of the boundary, in addition to generating an 
alternating term in S, have modified the uniform part of the entanglement entropy 
with respect to the result for periodic boundary conditions, Eq. (12. 2\\ . Note that both 
Eqs. (12. 3p and Eq. (12. 2\i only are valid for critical models. 

The impurity models we consider are defined using systems with open boundary 
conditions and when considering the impurity contribution to the entanglement care 
has to be taken with respect to the alternating term generated by the open boundary 
conditions. We do this by initially focusing only on the uniform part of the entanglement 
entropy and the contribution arising from the impurity to this part. Returning to our 
fundamental definition of Si mp , Eq. (II. 2p . we view S^no impurity) as the entanglement 
in a system without the impurity site, i.e. with one less site, R — 1 and all couplings 
equal to unity. We do not define S^no impurity) as S with J' K = since, as we shall 
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discuss in detail later, for J' K = the impurity spin can have a non-trivial entanglement 
with the rest of the chain. It then follows that the complete impurity contribution to 
the entanglement entropy cannot be obtained by subtracting results for a system with 
J' K = 0. We thus define the uniform part of the impurity entanglement entropy precisely 
as: 

S imp {J' K ,r,R) = S u (J' K ,r,R)-S u {l,r-l,R-l), r > 1. (2.4) 

For the subtracted part all couplings have unit strength. In a completely analogous 
manner one can also define the alternating part of the impurity entanglement entropy, 
Sfrnp, which we shall discuss in section [8] where it is shown that for the scaling 
form is modified from that of Eq. (11.61) . Our focus is therefore on the uniform part. 
For our numerical DMRG results a procedure for extracting the uniform and staggered 
part of the entanglement entropy is needed. We have found it sufficient to extract these 
functions by assuming that the uniform and staggered parts locally can be fitted by 
polynomials. If 7 sites surrounding the point of interest are used for the fitting a 7 
point formula can easily be derived as outlined in Appendix C] The numerical work 
has been done using density matrix renormalization group [22] (DMRG) techniques in 
a fully parallelized version, keeping m = 256 states. When performing calculations with 
J' K = and R even we have found it necessary to use spin-inversion symmetry [93] 
under the DMRG iterations in order to select the desired singlet ground-state. 

Other measures of the impurity entanglement entropy could have been defined and 
in section H] we will discuss one of these, the single site entanglement of the impurity 
spin with the rest of the chain, Si mp . As we shall see, this quantity gives an incomplete 
picture of the impurity entanglement. Another possibility would be to define a relative 
entropy [16] between the state with the impurity and a reference state without the 
impurity. It would be interesting to explore this latter possibility. We expect that 
similar scaling would be found as we show here with our definition of S imp . 

To gain some insight into how the impurity influence the entanglement entropy and 
lead to a non-zero S imp we show in Fig. [3] data for J' K = 0.41 for both R even, Fig. [3]^a) 
and R odd, Fig. [3](b). There are significant differences between the results for R even 
and R odd. It is also clear that the influence of the impurity is felt more strongly for 
R even compared to R odd. The resulting Si mp is clearly bigger for R even for all r 
as shown by the comparison in Fig. [3^b). In the limit where J' K — * this difference 
becomes more pronounced since Si mp for R even increases with decreasing J' K while for 
R odd it tends to zero with decreasing J' K . For the value of J' K = 0.41 used in Fig. [3] 
we shall later find that £k — 25.65 significantly smaller than R and it would have been 
natural to expect features in S imp signaling the presence of this length scale. However, 
Si mp is a monotonically decreasing function of r for both parities of R and no particular 
features are observed in S imp for r of the order of It is important to note that this 
fact does not imply a violation of scaling of the form Eq. (11.61) . 




Figure 3. (a) DMRG results for the total entanglement entropy, S(J' K ,r,R) for 
a 102 site spin chain at J|, with a J' K = 0.41 Kondo impurity (■) along with 
S(J' K = l,r — 1, R — 1) (□). Uniform parts (solid lines) and the resulting Si, np (J' K — 
0.41, r,R = 102) for R even, (b) DMRG results for the total entanglement entropy, 
S(J' K ,r,R) for a 101 site spin chain at J|, with a J' K = 0.41 Kondo impurity 
(■) along with S(J' K = l,r — 1,R — 1) (□). Uniform parts (solid lines) and the 
resulting Si mp (J' K = 0.41, r, i? = 101) for R odd. For comparison we also show 
Simp(Jfc = 0.41, r,R= 102) for R even from panel (a) (dashed line). 



J'k 


0.8 0.6 0.525 0.45 0.41 0.37 0.30 0.25 0.225 0.20 


£k R even 


1.89 5.58 9.30 17.40 25.65 40.5 111 299 ~565 ~1196 


£k R odd 


1.65 5.45 9.30 17.40 25.65 39.2 127 411 ~ 870 ~ 2200 



Table 1. The numerically determined values for £,k(J'k) using naive rescaling of 
Si m p{J'Ki r iR) a t nx ed r/R for J 2 = Jf. For R odd system sizes of R = 19 . . . 101 
have been used and for R even R = 18 . . . 102. The estimates become unreliable once 

6f » R- 



2.1. Scaling of Si mp at 

The scaling of Si mp (J' K ,r, R) for fixed r/R at J 2 = Jf was considered in detail in 
Ref. [75]. It was shown that Si mp with r/R fixed follows the expected scaling form, 
Eq. (11.61) . and is a function of a single variable r/^K- We expect this scaling to hold for 
r,R — r, £k ^> 1 when the results are not influenced by microscopic parameters such as 
the lattice spacing. For the moderate system sizes of R < 102 (R even) and R < 101 
(R odd) Simp shows a strong dependence on the parity of R and two different scaling 
functions for R even and odd were found. For r/^K ^ 1 these two scaling functions 
are essentially the same but they differ significantly for r/(K C 1. In the limit where 
R — > oo the two scaling functions eventually coincide. By requiring the data to collapse 
according to the scaling form, a naive estimate of can be obtained for both R even 
and odd. In tabled] we list the resulting £x obtained through such an analysis. 

We expect S imp to be a scaling function of only 2 parameters S imp (r /£ K ,r / R), or, 
equivalently Si mp (R/^K,r / R) for r,R,£,K ^> 1. It should therefore also be possible to 
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test this scaling by looking directly at S imp for fixed R/£k which should be a function of 
the single variable r/R. We use the previously determined £k listed in table[T]to test this 
assumption by selecting 2 sets of data, R = 400, J' K = 0.41, = 25.65, R/£k = 15.59 
and R = 86, J' K = 0.60, £ K = 5.58, R/£ K = 15.41. Our results are shown in Fig. 1 
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Figure 4. DMRG results for R = 400, J' K = 0.41, 6f = 25.65 yielding R/£ K = 15.59 
(open circles) and R = 86, J' K = 0.60, £, K = 5.58 with the ratio i?/6f = 15.41 (full 
triangles) plotted versus r/R. m = 256 states was kept in the DMRG calculation. 

where we observe an excellent data collapse using the previously determined values of 
£k confirming the expected scaling form. Combined with the results presented in [75J, 
this provides rather strong numerical evidence for the presence of the large length scale 
£k in the impurity entanglement entropy and demonstrates the scaling picture arising 
from the universal aspects of Kondo physics. 

2.2. The Fixed Point Entanglement 

It is instructive to consider what might be termed the fixed point entanglement, the 
entanglement occurring at the two fixed points at weak coupling J' K = and at strong 
coupling J' K — 1. When R is odd and J' K = the impurity spin is free while the 
remainder of the system is in a singlet ground-state. Due to the subtraction in Eq. (12.41) 
Simp as well as the alternating part is then zero. Surprisingly, as was shown in [75] , 
SimpiJ'x = 0; r i R) i s non-zero when R is even. We now discuss this in detail. 

We consider the case R even and the limit J' K — > 0. That is, out of the 4 degenerate 
states at J' K = we focus only on the singlet state, which is uniquely picked out by 
this limiting procedure. This state has a non-zero entanglement of the impurity spin 
with the rest of the system. On the other hand, to calculate the impurity entanglement 
entropy, defined in Eq. (I2.4|) . we must subtract S for a chain of odd length, R—l, with 
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J' K = 1. This odd length system has a spin doublet ground state. There is no simple 
relationship between S(r, R, J' K = 0) and S(r — 1, R — 1, J' K = 1). 

Explicitly, we may write the spin singlet ground state for a chain of even length R 
in the form: 

|^>=(l/v / 2)[|T>U>-|l>l^>]. (2.5) 

Here the first arrow refers to the state of the impurity spin and the second, double arrow, 
to the total S z of all the other spins. We may write the pure state density matrix as: 

1 



p P = \ifj >< ip\ 



T><T I ® I ^><^ I + I 1><I I ® I it><it I 
- | T><! I ® I ^><ft I - I !><T I ® I ftxfy 



(2.6) 



Now consider doing the partial trace over the spins at sites j with r < j < R. This trace 
can be done on each term separately in Eq. (12.61) . It is convenient to define 4 operators 
on the Hilbert space of the R — 1 spins 2, 3, ... R: 

Pm,m' = T^b \m >< m'\. (2.7) 
Here m =ff or Jj (i.e. ±1/2). Tr# means tracing over the spins at sites j with r < j < R. 

Pn = Pu- ( 2 - 8 ) 

While pjj4 and are themselves reduced density matrices for the R — 1 site model, 
Pfyft is clearly not a reduced density matrix since it is not Hermitian. The two density 
matrices p^ and p^ are related by spin-inversion. The total reduced density matrix, 
including the impurity spin, can then be written: 

p = -( Pu (2.9) 

2 \ ~Pu Pn J 

When considering the fixed point impurity entanglement entropy Si mp (J' k = 0,r,R) 
the term we subtract in Eq. (12.4j) . the entanglement entropy for a chain of odd length 
R — 1 with J' K — 1, is given by -Trp^ log p^ (or -Trp u log p w ). Since, is non- 
zero in the case at hand we see that Si mp (J' k = 0,r,R) is non-zero for R even, also in 
the limit R — ► oo. The resulting Si mp (J' k = 0,r/R) was numerically calculated at J| 
in [75] using DMRG methods and was shown to crossover in an approximately linear 
manner from ln(2) at small r/R to zero at r/R = 1. We show additional results in 
section [7J This cross-over can be understood [75] in terms of a term corresponding to 
the zero temperature thermodynamic impurity entropy In g shown to be present in the 
entanglement entropy at conformally invariant boundary fixed points [26] . 

If we now consider S imp (J' K = l,r/R) it will, up to a sign change, be the same for 
both R even and R odd since it is the difference in the uniform part of the entanglement 
entropy for an odd and even sized system. Since we expect the uniform part of the 
entanglement entropy to be independent of the parity of R in the thermodynamic limit 
we expect Si mp (J' K = l,r/R) to approach zero for increasing R as was demonstrated 
numerically in [75]. We show additional results for S imp (J' K = l,r/R) at J| in section [7J 
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3. Impurity Valence Bonds (IVB) and Single Particle Entanglement (SPE) 

In this section we define two heuristic quantities, the impurity valence bond (IVB) 
and the single particle entanglement (SPE). As we shall see these quantities capture 
essential parts of the impurity entanglement as described in the previous two sections. 
The IVB was first discussed in [75] and closely related ideas were developed by Refael 
and Moore in [91] (see also Ref. [95]). Here we recapitulate some of the results and 
further develop the ideas. In section [6] a much more complete development valid in the 
dimerized phase, J2 > Jf > showing more rigorously the presence of IVB and SPE terms 
in the entanglement entropy. 

We start by discussing the single particle entanglement. We think in terms of a 
tight binding model describing a finite chain with a single particle present in a state 
where the particle has probability p for being in region A and (1 — p) for being in 
region B. The wave-function can quite generally be written = £\ with \i) the 
coordinate space states, from which it follows that p = J2ieA l^| 2 - The reduced density 
matrix for region A can then be written p|l)(l| + (1 —p)\0) (0| where |1) is the state with 
the particle in region A and |0) the state with the particle in region B (absent from A). 
It immediately follows that the entanglement entropy is given by: 



We shall refer to this as the single particle entanglement contribution to the 
entanglement and we imagine that if a free spinon (or soliton) is present in the ground- 
state it should give rise to such a contribution to the uniform part of the entanglement 
entropy. Such a term would be present in the uniform part of the entanglement entropy 
for a uniform chain (J' K = 1) for R odd, where a single unpaired spin is present, but 
would likely be negligible for R even at J|. If we assume that the single spinon (soliton) 
picture is relevant at J| it follows that S im p(J' k = l,r/R) for R odd should be given by 
the SPE. In [75] this was shown not to be the case but as we shall see it the SPE is a very 
good approximation in the dimerized phase and in particular at the Majumdar- Ghosh 
point. 

The other component of the heuristic picture is the formation of an "impurity 
valence bond" (IVB) between the impurity and a site in the chain. See Fig. [5j When 
the IVB is cut by the boundary between regions A and B we expect it to give rise to a 
contribution of In 2 to the impurity part of the entanglement entropy. If the IVB does 
not cut this boundary the contribution to Si mp arising from the IVB is zero. We then 
see that: 



where p is the probability that the IVB connects the impurity spin to a site in region 
A. In the limit where J' k — > 0, the probability will, as above, be given by p = Yli&A l^i| 2 
where ipi now describe the wave-function of a single unpaired spin (soliton). At J£ we 
expect the fixed point entanglement entropy Si mp (J' k = 0, r/R) for R even to follow this 
form. In [75] this was shown to work well even at J2 = J 2- Care has to be given to 



Sspe = -plnp -(1-p) ln(l -p). 



(3.1) 




(1-p) In 2, 



(3.2) 
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the fact that the parity of R will influence this picture and we now discuss the different 
situations in a preliminary manner, postponing a more complete discussion to section [6j 




Figure 5. Typical impurity valence bond configurations for R even (a) and R odd 
(b). Note the unpaired spin on the 7 th site in (b). 



R even and J' K < 1: This situation is shown in Fig. 0(a). If J' K = we would expect 
the probability of forming an IVB between the impurity and a given site in the chain 
to be almost uniform throughout the chain. Hence, S imp should decrease approximately 
linearly with r/R from In 2 to zero as mentioned above. For J' K small but non-zero we 
expect the typical size of an IVB to be of the order of when £g- < R. In that case 
S irnp should decrease monotonically to zero for r/£ K ^> 1. For r /^ K) r/i?<C 1 we would 
expect it to approach In 2. This behavior is largely confirmed by the numerical DMRG 
calculations. 

R odd and J' K < 1: This situation is shown in Fig. [5](b). Since R is odd and the 
total spin of the system is 1/2 there is an unpaired spin present in the ground-state. 
When ^> R the unpaired spin is the impurity spin and we simply have Si rnp = 0. 
As £k decreases the probability of creating an IVB increases due to screening of the 
impurity. However, the average length of the IVB when it is present decreases with 
In section 12.11 we showed clear numerical evidence for a monotonically decreasing Si mp 
with r for R/£k constant, however, if we instead now imagine keeping r/R fixed and 
varying R/^k we see that the two above effects will trade off and give rise to a maximum 
in S imp for £ K « R [75] with R odd. 

Eventually, in the limit R — > oo, the parity of R no longer plays a role and we obtain 
the same S imp (r f° r both R even and R odd. However, for mesoscopic systems these 
effects could be of importance. 



3.1. The Entanglement entropy at the MG point 

While the entanglement entropy at the critical point, J|, is less amenable to the 
heuristic approach of this section, this intuitive picture sheds considerable light on the 
entanglement in the dimerized phase occurring for J 2 > J| as we now discuss. In 
fact, almost exact expressions for the entanglement entropy can be obtained for the 
Majumdar Ghosh model. In section [6] we present a much more detailed approach based 
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on a variational wave- function. At the MG point, the heuristic expression developed in 
this section for the entanglement entropy, follow directly from the variational approach. 



B 



(a) 



(b) 



Figure 6. MG state divided in two regions A and B. (a) The total number of site R 
is even and the ground-state is a singlet, (b) When R is odd, the ground-state contains 
a soliton f. 

We now focus on the Majumdar- Ghosh model, J2 = J/2. Consider first the case 
J' K — 1 and R even. Then the ground state is a trivial nearest neighbor valence bond 
state on every link between sites 2i — l and 2i, for i — 1, 2, 3 . . . R/2, as depicted in Fig. [6] 
(a). When r is even, S = since the ground state is a direct product of singlet states 
in regions A and B. On the other hand, for r odd, there is one valence bond, between 
sites r and r + 1, connecting regions A and B. In this case S = In 2, as exemplified in 
Fig. 0(a). 

Now consider J' K — 1 but R odd. The ground state now contains one soliton 
(unpaired spin) on one of the odd sites. This soliton can be in region A (see Fig. [6] (b)) 
with probability p(r) or region B with probability 1 — p, where: 

and i/jn° l describe the amplitude for the soliton to be on site n. The fact that regions 
A and B are sharing the soliton produces a "single particle entanglement" precisely as 
outlined above: 

S S pe= -p\np- (1 -p)m(l -p). (3.4) 

However, this is not the whole story. We must also take into account that there may 
be a nearest neighbor valence bond connecting regions A and B. Whether or not this 
is present depends both on the parity of r and on whether the soliton is in region A or 
B. If r is even, then this valence bond is present when the soliton is in A. Conversely, 
if r is odd, then it is present when the soliton is in region B. When this valence bond 
is present, it contributes an additional In 2 to S. The probability of it being present is 
p for r even and 1 — p for r odd. Adding this extra term we obtain: 

S(r) = -plnp- (1 -p)ln(l - p) + [1/2 + (-l) r (p - 1/2)] In 2. (3.5) 

As remarked earlier, when J' K = and R is odd, the impurity site is unentangled 
with the rest of the chain which contains only nearest neighbor valence bonds, between 
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sites 2i and 2i + 1. The only source of entanglement between A and B is a valence bond 
from site r and r + 1 when r is even. Thus 

S= (l/2)[l + (-l) r ]ln2. (3.6) 

Now, consider the case J' K = and i? even. Then there is an impurity valence bond 
stretching from site 1 to some even site. All other valence bonds have length 1. The right 
hand member of the impurity valence bond is again a soliton separating the 2 different 
nearest neighbor valence bond ground states, but it now forms a singlet with the spin at 
site 1. The soliton again contributes its single particle entanglement. Since the soliton 
forms a singlet with site 1 this contributes In 2 to S when the IVB terminates in region 
B, with probability [1 — p(r)] but make no contribution when it terminates in region 
A. Furthermore, there may be an additional nearest neighbor valence bond entangling 
regions A and B and contributing another In 2 to S. When r is even, this occurs when 
the IVB terminates in region B, with probability (1—p). When r is odd it occurs when 
the IVB terminates in region A, with probability p. Combining the contributions of the 
SPE, the IVB and the possible nearest neighbor valence bond gives: 

S = -plnp - (1 - p) ln(l -p) + [3/2 -p+ (— l) r (l/2 - p)\ In 2. (3.7) 

Both Eqs (13.71) and (13. 5p are shown to agree very well with numerical results at the 
MG-point presented in section [7J 



4. s imp and Weak Scaling Violations 

The simplest measure of how a qubit is entangled with the environment would be to take 
system A to be the impurity spin (qubit) itself and regard the environment as region 
B [53]. Often this is referred to as single site entanglement. If system A only contains the 
single impurity spin at the boundary of the chain it becomes impossible to extract the 
uniform part and we can therefore not use Eq. (12 Ah to define the impurity entanglement. 
Instead we have to use the complete entanglement entropy including both uniform and 
alternating parts. To distinguish it from Si mp in Eq. (12. we therefore denote it by Si mp 
and define it as S(J' K , 1, R) with no subtraction. We now analyze this quantity. 

First consider the case of R even. Then the ground state is a spin singlet. Since 
the impurity has spin-1/2, as does the rest of the system, we can write the spin- zero 
ground state in the form: 

|V >= (1/V2)[| T> ®| ^> -I l> ®| f>], (4.1) 

where the single arrow labels the state of the impurity and the double arrow the state 
of the rest of the system. Tracing out the rest of the system gives a two-dimensional 
density matrix for the impurity spin which is diagonal with elements 1/2 and hence a 
maximal entanglement entropy of In 2. The case of odd R is more interesting. Now 
the ground state is a doublet with total spin St = 1/2. Let us focus on the state with 
S T = +1/2. This must have the form: 

|V >= a\ T> ®|0,0 > +b\ |> <8>|0,1 >, (4.2) 
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where |0, > denotes a spin singlet state of the rest of the system and |1, 1 > denotes 
an S = 1, S z = 1 state of the rest of the system. All states are normalized to 1, so it 
follows that |a| 2 + |6| 2 = l. The density matrix is again diagonal with matrix elements 
| a, | 2 and \b\ 2 and hence 

s imp = -\a\ 2 \n\a\ 2 -\b\ 2 \n\b\ 2 . (4.3) 
On the other hand, the magnetization of the impurity in the ground states is given by: 

m imp = (H 2 - \b\ 2 )/2 = (2\a\ 2 - l)/2. (4.4) 
Thus we may write: 

)ln[(l/2±m imp )]. (4.5) 

± 

(This formula is also trivially true for the case R even in which case mj mp = 0.) For R 
odd, m imp , and hence s imp , shows an interesting dependence on R which reflects Kondo 
physics. m irnp was studied, for the usual fermion Kondo model, in [96] and (93, [98] for 
example. For weak Kondo coupling and relatively short chains, R -C mi mp ~ 1/2. 
On the other hand, for stronger Kondo coupling or larger chains, the magnetization 
is progressively transferred from the impurity to the rest of the chain, associated with 
screening of the impurity. In the limit of an infinite bare Kondo coupling, m imp = 
since the impurity spin then forms a singlet with one other electron and all of the 
magnetization resides in the other electrons, which have individual magnetization of 
order 1/R. 

Nonetheless, m imp is not a scaling function of R/£k- This is associated with the 
fact that the operator Si mp has an anomalous dimension [97], [98], ■jimp^Xx)- Thus it 
obeys the renormalization group equation: 



d d 
R ~dR + ^ K '~dX^ + limp ^ K ^ 



m imp (R, X K ) = 0. (4.6) 



The weak coupling /3-function, the variation of the effective Kondo coupling as we vary 
the length scale is: 

dXl< ~ -PM = X 2 k -Ix 3 k + ... (4.7) 



alliiL ' y ' "2 
We then see that if 7i mp (Ax) were zero, (14.61) would imply that function of 

X e ff(R) only, or equivalently a function of £,k/R only. However, a non-zero ■ji mp ^ 0, as 
occurs here, implies scaling violations. At lowest non-vanishing order [991193 EH] in A^, 

limp(\ K ) = Xj</2 + ... (4.8) 
The general solution of (14. 6 p is: 



m imp = exp{ / [-y(X K )/P(X K )]dX K }f(CK/R), (4-9) 
J o 

where / is some function of £,k/R or equivalently of Xk(R). From (14.71) and (14. 8 p we 
see that 

exp{ / [ 7 (A K )//3(A K )]dA K } = 1 + A° A -/2 + 0[(X° K ) 2 }. (4.10) 
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The fact that mj mp has this residual dependence on the bare coupling which cannot 
be adsorbed into the renormalized coupling at scale R implies a violation of scaling. 
However, since this effect vanishes as the bare coupling, — > 0, we refer to it as a 
"weak scaling violation" . 
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Figure 7. (a) Weak scaling violations for the single site impurity entanglement 
entropy Si mp (a) [Eq. (|4.5[) ] and the local impurity magnetization m imp = (Sf) (b). In 
panel (c) is shown an improved scaling plot of mi mp /(l + \® K /2). The £,k used in all 
panels are from table [TJ All data are for odd length chains between R = 19 . . . 101. 



This form of mj mp , exhibiting weak scaling violation, can be confirmed by an explicit 
perturbative calculation. This calculation was actually presented earlier in [96] for the 
ID lattice version of the fermionic Kondo Hamiltonian: 

R-2 

H = -t^(^% +1 + h.c.) + J K S imp 

3=1 

In [96] R — 1 was called L. The result was: 

2 



t a 



(4.11) 
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(L + 1)J 
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k,k' 



0{e k ,)0{-e k ) 



sin k sin k' 



(4.12) 



Here 6(x) is the Heavyside or step function, and e& = — It cos k is the dispersion relation. 
The allowed values of k and k' occurring in the sum in f)4.12p are k = nn/R for 
n = 1,2,3,.... Assuming R ^> 1, we replace the sums by integrals. The integrals 
diverge logarithmically at k = k' = kp = ir/2. Letting k = ir/2 — q and k' = ir /2 
the integrals near q, q' — take the form: 
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1 
2 



J 



~1 2 



7Tf 



ln(i?/a). 



(4.13) 



Here we have replaced the small g limit of the integral by 1 /R and the large q limit by 
a short distance cut off, a, of order a lattice spacing. The continuum limit of this tight 
binding model gives the Hamiltonian of (1A.3[) with Ax = Jx/^t and v = 2t, so we may 
write this as: 

2 

\n(R/a). (4.14) 



m 



imp 



The lowest order correction to the effective coupling as determined by the first 
(quadratic) term in the /^-function is: 

X K (R) ^\° K + (\° K ) 2 HR/a 1 ) + ..., (4.15) 

where a' is another short distance cut off. Thus we see that, to 0[(A^) 2 ] we can write: 
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imp 
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This has the form of (14. 9p with the scaling function: 
flX K (R)] * \ - • ~ 1 

The . 



(4.16) 



(4.17) 



4 ••"2 41n(i?/60 
represents terms of 0(Xk(R) 2 and higher in the effective coupling at scale R. 



The fact that m; mp exhibits weak scaling violations implies that Si mp does also, as 
displayed in Figs. [3(a). In Fig. [7] we show DMRG results obtained at the critical point 
J 2 = 0.2412 for various odd lengths 19 < R < 101 and Kondo couplings 0.2 < J' K < 0.8. 
The local magnetization at the impurity site m imp = (S^) is shown in Fig. [3(b) as a 
function of R/£k, using the values of £k(J'k) determined previously (tabled]). Clearly 
both Si m p and mi mp violate scaling since the various curves do not fall on top of each 
other. As outlined above we expect m imp /(l + X° K /2) to scale much better and this is 
shown in Fig [7(c). Some deviations from scaling are still visible but clearly the scaling 
has improved. We expect the remaining discrepancies could be improved upon by using 
more optimal values for £x instead of the ones determined from other scaling plots which 
have significant uncertainties associated with them for either very large or very small 

The Kondo physics also is expected for J 2 < J| and therefore also for the 
unfrustrated chain at J 2 = [55]. The unfrustrated spin chain model can be 
investigated using Quantum Monte Carlo (QMC) methods as an alternative to the 
DMRG computations. Since QMC works for any temperature, it allows to investigate 
finite temperature scaling properties like the spin susceptibility (see Ref. [55]) from where 
we can get precise estimates for the Kondo temperature T K = v s /£k- These estimates, 
given in table [2] (along with estimates obtained in Ref. [55] by solving the Bethe Ansatz 
equations for this model |100j ). are then used to check the scaling violations at J 2 = 
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J'k 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.8 




~ 4 x 10 3 


168 


37 


14 


6.24 


3.5 


1.8 


0.95 



Table 2. £k(J'k) for J 2 = 0, as determined from scaling of the finite temperature 
impurity susceptibility and Bethe Ansatz calculations (reported in Rcf. [55]). 



of the local magnetization at the impurity site as well as the single site impurity 
entanglement entropy. As in the frustrated case discussed above, also here the weak 
scaling violations are clearly present. In Fig. [HI we show the QMC results obtained for 
various odd lengths 9 < R < 257 and Kondo couplings 0.1 < J' K < 0.8 where mi mp 
violates scaling as well as s imp . 
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Figure 8. Weak scaling violations for the local impurity magnetization mj TOp = (5f ) 
(inset) and the single site impurity entanglement entropy Si rnp [Eq. (|4.5|) ] (main panel). 
QMC results obtained for the critical spin chain [model (|1.4|) ] at J2 — with various 
odd lengths R and Kondo couplings J' K indicated on the plot. The scaling violations 
are shown as a clear non-universality as a function of R/£k(J'k) 



4-1. Discussion 

In [971 [98] a rather general discussion was given of which physical quantities are given 
by pure scaling functions and which ones exhibit weak scaling violations. In general 
quantities defined far from the impurity, such as the local susceptibility, are pure scaling 
functions whereas those that depend explicitly on the impurity spin operator exhibit 
weak scaling violations. In particular, the impurity susceptibility and specific heat or 
(thermodynamic) entropy are given by pure scaling functions. This is because they 
only involve conserved quantities, the total spin and Hamiltonian operators which have 
zero anomalous dimension. (Roughly speaking, for the case of the total spin operator, 
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the anomalous dimension of the impurity spin operator is canceled by the anomalous 
dimension of the electron spin density operator near the impurity.) We expect the 
impurity entanglement entropy, Si mp (J' K , r, R) to be a pure scaling function when r 1, 
since it is a quantity defined far from the impurity. Indeed, in the approach of Cardy 
and Calabrese, reviewed in Sec. [5] the entanglement entropy is obtained from Green's 
functions of an operator $ n , which is inserted at the positions ±r. Since this operator 
is inserted at a location far from the boundary, we expect that it will not have any 
anomalous dimension associated with the Kondo interaction. Therefore the one-point 
function obeys the RG equation of Eq. (I4.6p with zero anomalous dimension, implying 
that it is a pure scaling function. Indeed our DMRG data is consistent with such scaling 
behavior as we have shown in Fig. (j3J). 



5. Fermi Liquid Theory for S,, 



imp 



In the limit £x <C r, within the FLT as outlined in Appendix A, we can also 
calculate Si mp by treating H int in Eq. (1A.17I) in lowest order perturbation theory. The 
entanglement entropy is obtained by means of the replica trick. If Trp(r) n is known, 

S= -lim-^-[Trp(r) n ], (5.1) 

n—>i an 

where p(r) is the reduced density matrix for the subsystem (0, r). In the path integral 
representation of Euclidean space-time, 

Trp(r)" = (5.2) 

where Z n (r) is the partition function on an n-sheeted Riemann surface TZ n , with the 
sheets joined at the cut extending from r to R [21 (2Sj • Now the original problem has been 
transformed into the calculation of the partition function with a nontrivial geometry. 
We use the approach where the Hamiltonian is written in terms of left movers only, 



obeying PBC on an interval of length 2R. (See Appendix A). In the critical region, the 
system is conformally invariant and CFT methods are applicable. Starting with zero 
temperature, the n-sheeted Riemann surface 7Z n , can be mapped to the usual complex 
plane C [26j . Then the expectation value of the energy momentum tensor T on TZ n is 
simply given by the Schwartzian derivative: 

jrpf y, A n ( Ul -U 2 ) 2 

(T(w)) nn = — -, (5.3) 

where A n = (c/24)(l — (1/n) 2 ) and we have u\ = ir and u 2 = —ir here. Then Calabrese 
and Cardy [26] observed its important connection to the correlators on C through the 
Ward identity: 

J[d(f)}T(w)e~ SE ^ 



{T(w)) 



~ f[d(t)]e- s v(Kn) 

(5.4) 



(T(«;)$ B (r)*_ B (-r)) c 



($ n (r)$_ n (-r)) f 
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The fictitious primary operators $± n on the branch points have the left scaling 
dimensions A n . They concluded that Trp n = Z n /Z n behave identically to the n-th 
power of ($ n (r)$_ n (— r)) c under the conformal mappings or explicitly, 

Trp n = £ n (2r/a) (c/12)(n " 1/n) . (5.5) 

Applying the replica trick, the entanglement entropy is S ~ (c/6)ln(2r). Then they 
extended the result to finite system size R or infinite system size and finite temperature 
T = 1/(3 by applying the corresponding conformal mapping to ($ n (r)$_ n (— r)) c |26j . 

Now with the presence of the local irrelevant interaction Eq. ( 1A.17I) . we should 
calculate perturbatively the correction to the partition function Z n in order to get the 
impurity entanglement entropy Si mp . Luckily, the irrelevant interaction is just the energy 
momentum tensor itself and its expectation value on the n-sheeted Riemann surface is 
just Eq. ( 15. 3p . The correction to Z n of first order in £ K is: 

/oo 
dT(H s , L (T,0))n n - (5.6) 
-oo 

Ti s ,L = T/(2ir) where T(r, x) is the conventionally normalized energy- momentum tensor 
for the c = 1, free boson conformal field theory corresponding to the spin excitations of 
the original free fermion model. After doing the simple integral and taking the replica 
limit, for R — ■> oo, we get 

5^ = 7^/(127-). (5.7) 

In principle, in order to extend the FLT calculation to finite R, we will need the 
conformal mapping from a infinite n-sheeted Riemann surface to a finite one. On the 
other hand, we can also try to exploit Eq. (15 ,4p following ideas similar to Ref. (26] by 
applying the standard finite size conformal mapping to both (T$ n $_ n ) c and ($ n $_ n ) c . 
Then, in the first order perturbation, 

- 6&) = -(-)- f dr[(^) p^ILim f , (5.8) 



We use the integral 



csct, (5.9) 



cosh ax — cos t a 

from Ref. |101j and differentiate Eq. (15.91) with respect to t on the both sides. Applying 
this result and the product-to-sum hyperbolic identity to Eq. (15. 8p . we can complete the 
integral and after the replica limit we get: 

S imp = ^[1 + tt(1 - -) cot(— )]. (5.10) 

Of course, Eq. (15. lOf) reduces to Eq. (15. TP for R ^> r and both of them agree with the 
scaling form of Si mp . Interestingly, Eq. (1 5 . 1 [) can be regarded as the first order Taylor 
expansion in £ K /r and £k/R of Sjj = (1/6) ln[i?sin Trr/R] with r and R both shifted by 
7r^A-/2. Consistently, Eq. (15.71) can also be obtained from expanding (1/6) ln(r + 7r^/2). 
In fact, many other quantities such as impurity susceptibility, specific heat and ground 
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state energy correction can be also obtained in this fashion by shifting the size of the 
total system, R, to R + irt; K /2. 

Within CFT methods, we can also calculate Si mp for infinite R but at finite 
temperature (3. We apply the standard finite temperature conformal mapping to 
(T$ n $_ n ) c and ($ n $_„) c . The result for first order perturbation is just to replace 
2R by (3 and sinh by sin in Eq. (15.81) with the integral from —(3/2 to (3/2. Completing 
the straightforward integral yields 

S imp = {n 2 Z K T/(6v)] coth(27rr:Z», (5.11) 

valid for T,v/r -C Tk- In the intermediate temperature regime, v/r <C T <C Tk 
and hence rT ^> v, Eq. (15. lip approaches the thermodynamic impurity entropy, 



S, 



imp 



tt 2 UT/(6v) = n 2 T/(6T K ) 



Cimp\ 



(T), the well-known impurity specific heat. 
(Tk = v/£k is the Kondo temperature.) This is consistent with the observation that 
in this limit the entanglement entropy approaches the thermodynamic entropy as noted 



in Ref. [26]. In Appendix E this connection is explored in more detail and it is argued 
that quite generally S will approach the thermodynamic entropy for T ^ v/r. 



[1+7[(1-l7R)cot(7tl7R)]7i/12 

S imp < J 'K= 1 - ' r . R = 40 °UK= - 64 

□ S (J' K =1.0,r,R=401U K =0.65 




400 



Figure 9. Si„ lp (J' K = l,r,R) for R = 400,401 compared to the FLT prediction, 
Eq. (I57L01) . 

The FLT expression for the uniform part of the entanglement entropy, Eq. (15.101) 
has been compared in detail with numerical DMRG data in Ref. [75]. The resulting 
values for are shown in table [31 Excellent agreement between this expression and 
numerical data for Si mp (J' K ,r, R) was found for J' K < 1. However, the fixed point 
impurity entanglement entropy Si mp (J' K = l,r,R) with J' K — 1 itself should in fact be 
given by this same expression for R both even and odd. Now, with a = 0.65 of order 
(1). In Fig. Owe show results for S imp (J' K = l,r,R) for R = 400,401. In both cases 
do we find good agreement with the FLT result, Eq. (15.101) . The small discrepancies 
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J'k 


1.00 0.80 0.60 0.525 0.45 0.41 0.37 0.30 




0.65 1.97 5.93 9.84 17.83 25.65 38.29 83.79 



Table 3. £k(J'k) as determined from S imp (J' K ,r, R) for R = 400 using the FLT 
prediction, Eq. (|5.10p . The estimates become unreliable once becomes comparable 
to R. 

between the numerical data and the FLT results for r close to 400 are likely due to 
the approximate forms used for extracting the uniform and alternating parts of the 
numerical data. For J' K ^ 1 these discrepancies are less pronounced since the numerical 
signal for Si mp is larger by a factor of R/£k- 

It is important to note that the observed agreement implies that the fixed point 
impurity entanglement entropy, S imp (J' K — 1, r, R) ~ 1/R for both even and odd R and 
hence is zero in the thermodynamic limit. 

6. The TS- and OTS-ansatz for the Dimerized Phase, J 2 C < J 2 < J/2 

For J2 > J 2 the spin chain model, Eq. fll.4j> . develops a gap and a non-zero dimerization. 
The well known Majumdar- Ghosh (MG) model [58] (J2 = J/2) is part of this phase and 
we therefore take the MG model as our starting point for an analysis of the entanglement 
entropy in the dimerized phase. The MG point constitutes a disorder point [102] beyond 
which the short-range correlations become incommensurate and our analysis therefore 
focus on the regime J 2 < J2 < J/2. Due to the gap above the dimerized ground- 
state, it becomes possible to proceed using variational wavefunctions for this range of 
parameters, yielding surprisingly precise results. The initial assumption is that the 
lowest excitation in the dimerized phase is described by a single free spin acting as 
domain wall or a soliton. This is an approximation since a complete description would 
include states where the domain wall extends over several sites. The thin soliton states 
are not orthonormal. By assuming that they are, we arrive at a further simplified ansatz 
which never the less proves to be surprisingly precise. Below we detail this approach. 

6.1. The Thin Soliton (TS) Ansatz 

For the MG model [58J (J2 = J/2) h is well known that for R even the singlet ground- 
state is two fold degenerate [58], [91] I103[ [104J. These two ground-states corresponds to 
the formation of singlet states between nearest neighbor spins either between sites 2n + 1 
and 2n + 2 or 2n + 2 and 2n + 3. The resulting dimerization then occurs in two distinct 
patterns. (See Fig. [6]) For R odd a very good approximation to the ground-state is 
obtained by assuming that a single spin, the soliton, is left unpaired separating regions 
with dimerizations in the two above mentioned patterns. If we number the sites of the 
system r = 1, . . . , R the number of odd sites is given by N° = (R+ l)/2. We define a 
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state \n) by: 

|n) = |^".~^ T -•••-). (6.1 

Here, — indicates a singlet between site r and r + 1 with n = 0, . . . = N° — 1 
such singlets occurring before the soliton indicated by the f- Here, n can take on 
the values n = 0, . . . = N° — 1 with the total number of dimers. Note that, 
singlets to the left of the soliton occur between sites [2n + 1, 2n + 2] and between sites 
[2n + 2, 2n + 3] to the right of the soliton. We use the Marshall sign convention that 
[2n + 1, 2n + 2] = (| tl) - | IT))/V2 whereas [2n + 2, 2n + 3] = (| |t) - | U))/V2. 
We shall refer to these states as thin soliton states (TS-states) since the soliton is not 
"spread" out over several sites by including valence bonds of more than unit length. 
Note that, in such a dimerized state the soliton can, for R odd, only be situated on odd 
sites r = 1, 3, 5, ... R. We can then write an ansatz for the ground-state wavefunction 
in the following manner [9H HHl HUH [T05] : 

N d 

|^ 5 >^E^°V>- (6-2) 

n=0 

We shall refer to this ansatz as the thin-soliton ansatz (TS-ansatz). If we consider R = 5, 
we find 3 linearly independent (but not orthogonal) thin soliton states. However, with 
5 S = 1/2 spins it is easy to see that there are in fact 5 5=1/2 states. The thin 
soliton states do therefore not form a complete basis for the S = 1/2 subspace and the 
TS-ansatz is therefore variational in nature for R odd. However, for the MG-model 
it is known that the TS-ansatz is very precise [1031 11041 1105] . Caspers et al [1031 1104] 
improved on the TS-ansatz by including terms with longer valence bonds in a systematic 
manner and showed that the resulting variational energies were only changed slightly. 

For R even Shastry and Sutherland [91] considered excited states corresponding 
to 2-soliton states and studied bound-states of solitons with relatively high energies. 
Subsequently, exact wave-functions for bound soliton states were found by Caspers and 
Magnus [103] . However, at low energies the solitons behave as free massive particles |105] 
with the soliton mass defined by: 

A sol = lim E(R) - E ( R +V + E(R-V ^ R Qdd _ (g 3) 

R—fOo 2 

Here E(R) is the ground-state energy for a system of length R. Using DMRG the mass 
of the soliton has been estimated |105] for the MG model: 

^oi IRG /J = 0.1170(2). (6.4) 

Using the TS-ansatz and periodic boundary conditions (pbc) the soliton mass was 
determined to be [1031 HU1] A^' pbc /J = 0.125 and including 3 and 5 spin structures in 
the variational calculation [1031 1104] the estimate improved to 0.11701. Note that the 
ground-state energy, E(R), is an extensive quantity and the term cx R, 3JR/8, is given 
exactly by the TS-ansatz. The error in the estimate of the ground-state energy using 
the TS-ansatz is only (0.125 — 0.11701) J ~ 0.008J, a small quantity independent of R. 
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It is useful to obtain an estimate of the wave-function for a single soliton from a 
simple physical picture. Such an estimate can be obtained in the following manner: The 
soliton is repelled by the open ends used in the present study and we therefore expect 
the thin soliton to behave particle in a box [105] with ijj s _°l = = 0. In that 

case we find: 



C°^A/l^sin(^±^). (6-5) 



N d + 2 \ R + 3 

The ip'™ 1 can also be determined using variationally methods as shown in Appendix F 
and in the following we shall use such variationally determined . However, as shown 
in Appendix F the variational estimate is only marginally better than the above form, 
Eq. (16.51) . Several other quantities can also be calculated for the MG model using the TS- 
ansatz, such as the on-site magnetization, (S*), and the spin-spin correlation function, 
(S r ■ S r+ \). Since these quantities are unrelated to the main focus of the present paper, 
the entanglement entropy, we have included them in Appendix F 



We now turn to calculating the entanglement entropy for the Majumdar- Ghosh 
model using the TS-ansatz. The two cases of R even and odd have to be considered 
separately and we start with R odd. 

6.1.1. TS-ansatz for the Entanglement entropy for R odd, J' K = 1 Using the TS-ansatz 
it is possible to obtain an explicit expression for S(J' K = l,r,R) for R odd. We start 
by separating Eq. (16.21) into contributions from region A (1 . . . r) and B (r + 1 . . . R). 
Without loss of generality we initially assume that r is odd. We write: 

3 

|¥Js> = £ C ^>l&>> ( 6 - 6 ) 

i,j=0 

with the states in the Hilbert space for region A and the \(f)j) states in the Hilbert 
space for region B. We find for 

2 



n=0 



i — 1 
2 



iifc> = r ••• — ' r> 



i — i 

2 



llfe> = l' ••• ' !)• (6-7) 

For \<j)j) we then write: 
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= E t r 



n=0 



With these definitions we find: 



n 5 ) = l^l)l0l) + -y|[l^)l02)-^3)|03)] 



and thus: 



C 



( i o \ 

A 



-\ 



(6.8) 
(6.9) 

(6.10) 



V2 I 



Clearly, having defined the states this way, the and <pj) are not orthonormal. It is 
easy to see that (ip^lipi) — for i — 1,2 and trivially (^2^2) = (^3^3) = 1; however, 
using the Marshall sign convention described at the start of section 16.11 we find, 



r-l 
2 



-n\ 



n=0 

r-l 



■m| 



(6.11) 



n,m=0 



Likewise, we see that (0 3 |0j) = for j — 1, 2 and (0i|0i) = 1, however, 



'1^2 



n 7 _ ^ 

v n=0 



^2 VI/ 



-|n— m| 



n,m=0 

= (0 3 |0 3 ). (6.12) 

In order obtain the reduced density matrix, p, describing region A with the 
usual properties (i.e. Trp = 1 with eigenvalues uJi > 0) it is easiest to proceed by 
orthonormalizing the states This is done in the usual manner by defining: 

With this orthonormalization the coefficient matrix C then becomes: 

C= ^(^2|^2)-|(V2|^l)| 2 ) 





IV2)- l^iX^I A) 



(6.13) 



V 



(6.14) 
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The reduced density matrix for region A can now be obtained by tracing out region B: 

3 



jl ,32=1 



(6.15) 



from which the entanglement entropy is calculated employing the standard formula 
S = — Trplnp. An equivalent expression for r even can be obtained by interchanging 
\i>i) and \4>j). 

As an example of how well the TS-ansatz, Eq. (16. 2p . works we show in Fig. [TOlthe 
results of a calculation of S( J' K = l,r,R— 201) using this ansatz compared with DMRG 
results for the same quantity. The agreement is almost perfect. In Fig. [TOl variationally 
determined ^ ol ( Appendix F ) have been used. 



T — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I - 



W 0.6 



o .oS(r K =l,r,R=201) 
• TS-ansatz 




200 



Figure 10. DMRG results for S(J' K — l,r,R= 201) for the spin-chain model with 
J2 = J/2 (MG model). The small circles represent the theoretical result obtained 
using the ansatz Eq. (|6.2[) with variationally determined 4>n° l ( | Appendix F[ ). 



6.1.2. TS-ansatz for the Entanglement entropy for R even, J' K = It is now straight 
forward to extend our results from the previous section to the case where R is even 
and J' K = 0. Our starting point is the expression Eq. (12. 9p for the reduced density 
matrix. In order to calculate and in particular we need an ansatz for with 
all couplings uniform and total length R = R — 1 odd with r' — r — 1 denoting the 
spatial coordinate. This can straightforwardly be obtained from Eq. (16.21) by inverting 
the soliton spin to obtain: 

|$} ff )~^|^ I -...-)• (6.16) 

n=0 

Using this ansatz we can now evaluate as well as p^. 
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As above, we start by separating Eq. (I6.16P into contributions from region A 
(1 . . . r') and B (r' + 1 . . . B!) again initially assuming that r' is odd. We write: 

3 

Hs) = Y< C ^'M), (6-17) 

i,j=0 

with the |^) states in the Hilbert space for region A and the |<//) states in the Hilbert 
space for region B. In this case it is convenient to define |^): 



r'-l 
2 

oh 



n=0 



r'-l 
2 



i^) = r ... — * i) 



r'-l 
2 




i^) = r ... — ' t). (6.i8) 

For we write: 



102) 



n=0 

(6.19) 

With these definitions one finds that (ipi\ipj) = (ip'^ipj) and (4>i\4>j) = (<fiiWj)- Hence, 
for this choice of states, p^ is numerically identical to p-^, Eq. (I6.15p . However, the 
actual states \<p'j) are not identical to the and (<fii\<f>'j) is not the identity matrix. 
This becomes important for the evaluation of p-^jj. which can be expressed: 

3 

31,32=1 

with the coefficient matrix C given as before by Eq. (I6.14p . It is easily seen that 
(<t>iWj) = except for = (falfa), (fcWi) = (<fa\<f>i) and (<j>i\<f/ 2 ) = (<f>i\(fa). p U 

can then be evaluated along with p^ = pin and the full density matrix, Eq. (12. 9p . 
describing J' K = and B even constructed. Again, by essentially interchanging the role 
of |^) and \4>j), it is straightforward to obtain results for r' = r — 1 even as well. 

With the matrices, pff, p^, Pf^piyf in hand, the full density matrix for J' K = 
and B even can be constructed from Eq. (I2.9p and the complete entanglement entropy 
calculated. 

As an example of how well the TS-ansatz, Eq. (16. 2p and Eq. (I6.16p . work we show 
in Fig. [LTJ the results of a calculation of S(J' K = 0,r,B = 200) using this approach 
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Figure 11. DMRG results for S(J' K = 0,r,R= 200) for the spin-chain model with 
J2 = J/2 (MG model). The solid line represents the theoretical result obtained 
using the TS-ansatz Eq. (|6.2p and Eq. (|6 . with variationally determined ipn" 1 
QAppendix~T| . 



compared with DMRG results for the same quantity with variationally determined i\) s ° l 
(Appendix F). Excellent agreement is observed. The small discrepancies between the 
DMRG results and the TS-ansatz visible at a few r values are due to complications 
using spin-inversion in the DMRG calculations specific to this value of J 2 . 

We end this section by remarking that in principle the TS-ansatz could be used for 
any J 2 > J 2 by simply redoing the variational calculation determining the ipn 1 - 



6.2. The Orthonormal Thin Soliton (OTS) Ansatz 

While relatively straight forward to use, the TS-ansatz does not yield expressions that 
are intuitively easy to grasp. It therefore seems desirable to develop a simplified picture 
of the underlying physics which we will now try to do. We shall do this by assuming 
that the TS states are orthonormal, thereby arriving at a simplified ansatz. Within 
this simplified ansatz the "single particle entanglement" and "impurity valence bond" 
contributions to the entanglement entropy, introduced in in sub-section [31 are derived 
explicitly employing this ansatz. 



6.2.1. OTS-ansatz for the Entanglement entropy for R odd, J' K = 1 We begin by 
focusing on the case of R odd and J' K — 1 at the MG-point (J 2 = J/2). Our 
starting point for the calculation of this quantity using the TS-ansatz was the expression 
|\E^ >= jCij\tpi)\(t)j) . We now make the simplifying assumption that all the \ipi) 
and \(j)j) are orthonormal. We shall refer to this as the orthonormal thin soliton ansatz 
(OTS-ansatz). This does not change the coefficients Cy: C^i = 1, C 2>2 = —£3,3 = l/y/2 
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and otherwise. However, the reduced density matrix now significantly simplifies and 
assuming orthonormality we find for r odd: 

/£n=ol^T o \ 

Pn= lEn4K°f o . (6.2i) 

Here, k = (r + l)/2 and, as above, = (R — l)/2. If we by p denote the probability 
that the soliton is to the left of the point r we then see that this is simply: 

\ 



Pn 



o 

2 






1-p 

2 



r odd. 



(6.22) 



/ 



It immediately follows that for r odd: 

S(J' K = l,r,R) = -pln(p) ~ (I-P)ln(l ~P) + (l-p)ln(2). (6.23) 
If we now turn to r even while still considering R odd and = 1 we instead find: 
/ 1 -p \ 

(6.24) 



Pn 



\ 



o 
o 



2 

2 u 

g 



r even. 



This naturally follows from the fact that now describes a situation with the soliton 
to the right of the point r while p% 2 and p% 3 have the soliton the left. We then see that 
in this case for r even: 

S(J' K = l,r,R) = -p\n(p) - (l-p)ln(l -p) +pln(2). (6.25) 

We emphasize that Eq. (16.231) and (I6.25P precisely equalise heuristic expression Eq. (13.51) 
based on the SPE and IVB. 

We can now explicitly obtain expressions for the uniform and alternating part of 
the entanglement entropy, we find for R odd: 



Su(J' K = 1, r, R) = - p ln(p) -(1-p) ln(l ~ p) + \ ln(2) 
S A (J' K = l,r,R) = ( 1 --p)ln(2). 



(6.26) 
(6.27) 



From this we can immediately extract Si mp (J' K = l,r,R) for R odd since at the MG- 
point Su{J'k — 1; r i R) for R even is simply ln(2)/2 for any r. We then find: 

S imp {J' K = 1, r, fl) = -pln(p) - (1 - p) ln(l - p), i? odd. (6.28) 

Since this is effectively the entanglement resulting from the presence of a single thin 
soliton in the ground-state we see that the Si mp in this case is given uniquely by the 
SPE. 

It is of considerable interest to test how well the TS- and OTS-ansatz agree at 
the MG-point and to test how well either one of them agree with the DMRG results. 
In Fig. fl2Ta) we show results for the difference between S imp (J' K = l,r,R) calculated 
with DMRG and the TS-ansatz scaled by R. Fig. [T2T b) shows results for the difference 
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Figure 12. Comparison of results for Si mp (J' K = l,r,R) at the MG-point J2 = J/2 
for a range of system sizes obtained using DMRG, TS-ansatz and the OTS-ansatz. (a) 
Difference between DMRG results and TS-ansatz results scaled by R. (b) Difference 
between OTS and TS-ansatz results scaled by R. For the TS-ansatz variationally 
determined ( [Appendix F[ ) have been used where as for the OTS-ansatz the 

analytical result for p(r) Eq. (|6. 38[) have been used. 



between S imp (J' K = 1, r, R) calculated with the OTS-ansatz, Eq. (I6.28p . and with the TS- 
ansatz, again scaled by R. The data in Fig. [T2l clearly show that in the thermodynamic 
limit all three approaches agree. This somewhat surprising result indicates that at the 
MG-point the non-orthogonality of the TS-states for the calculation of the entanglement 
entropy cannot play an important role. 



6.2.2. OTS-ansatz for the Entanglement entropy for R even, J' K = We now turn to 
the case where J' K — and R is even. We then have to consider the full 6x6 density 
matrix. If we initially take r even we find (R even): 
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(6.29) 



This follows from the observation that within the OTS ansatz and with the states 
and \4>i) as defined in Eqs. (16.81) and ( 16.191) . (4>k\4>'i) = unless k — I — 1 in which case 
(<t>i\4>i) = 1- By diagonalizing this matrix it follows that for r even (R even): 

S(J' K = 0,r,R) = -pln(p) - (1 -p)ln(l - p) + 2(1 -p)ln(2). (6.30) 
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In the same manner we write for r odd (R even 
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(6.31) 
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In a way analogous to Eq. (16.291) . we find non-zero off-diagonal elements. From which 
it follows that for r odd (R even): 

S(J' K = 0,r,R) = -p\n(p) - (1 - p) ln(l - p) + ln(2). (6.32) 

Again we emphasize that Eq. (16.301) and (I6.32p precisely equal the heuristic expression 
Eq. (E2D based on the SPE and IVB. 

As above, we can now explicitly obtain expressions for the uniform and alternating 
part of the entanglement entropy, we find for R even: 



Su(J'k 
Sa(J' k 



0,r,R) 
0,r,R) 



pln(p) -(l-p) ln(l -p) + (--p) ln(2) (6.33) 



l^-P) ln ( 2 )- 



(6.34) 



We note that the uniform part is simply SPE +(3/2— p) ln(2). Again we can immediately 
extract Si mp (J' K = 0,r,R) for R even at the MG-point by using the above result for 
S V (J' K = l,r,R), Eq. fl6\26|) . We find: 

S imp ( J' K = 0, r, R) = (l- p) ln(2), R even. (6.35) 

We note that this expression does not contain a contribution from the single part 
entanglement (SPE), but purely a term related to the impurity valence bond (IVB). 

By numerically calculating the density matrix p for large R using the TS-ansatz, 
we have verified that the p indeed does follow the above approximate forms, Eqs. (16.291) 
and (16. 31 j) . for both r even and odd. 

It is useful to have an analytical expression for p, the probability of finding the 
soliton in region A (x < r). Simplifying our expression for the soliton wavefunction, 
Eq. (16.51) . we write: 
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(6.36) 



With p corresponding to the probability of finding the particle in region A (x < r) we 
can then write: 
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From which we find: 
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(6.37) 



(6.38) 
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This expression agrees rather well with the probability extracted from the variationally 
determined soliton wave-function. See Appendix F Since p only depends on the single 
variable r/R we see that also Si mp (J' K = 0,r,R) and Si mp (J' K = l,r,R) as obtained 
from the OTS-ansatz are functions of the single variable r/R. When we refer to the 
OTS-ansatz we always assume that p(r) has been determined using the above analytical 
form, Eq. (KT5g|) . 

If a similar analysis had been performed away from the MG-point, at a different 
J 2 > J%, the soliton wave function, Eq. (16. 36ft would likely change and thereby the 
above expression for p. However, it is noteworthy that the remaining expressions such 
as Eqs. (16. 28ft and (16.35j) only depend on J 2 through p. 



7. Numerical Results for the Fixed Point Entanglement 

7.1. Fixed Point Entanglement at the MG-point, J% = J/2 




Figure 13. DMRG results for S imp {J' K = l,r,R) for R = 101,401 for the spin-chain 
model with J2 = J/2. The dashed line represents the OTS-ansatz result Eq. (|6.28p . 
For R = 101, 401 we also show the results for the TS-ansatz Eq. (|6.2[) as the solid and 
dotted lines. 



At the MG-point (J2 = J/2) we can use the TS-ansatz, Eq. (16.21) . directly to 
calculate Si mp (J' K = l,r,R) for R odd. This is straight forward since S(J' K = l,r,R) 
for R even is simply ln(2) for r odd and zero otherwise. The result of such a calculation 
it the MG-point is shown in Fig. [12] (solid lines) for R = 101,401 where we also show 
result for the same quantity calculated using the TS- and OTS-ansatz, Eqs. (16.21) . with 
variationally determined ^ s ° l \ and (16.281) with p{r) from Eq. (I6.38p . Note that, the 
TS-ansatz depends on R where as the OTS-ansatz is independent of R. Excellent 
agreement is observed between the theoretical and numerical results. Both the DMRG 
results and TS-results rapidly approach the OTS result as implied by Fig. [12j Clearly 
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this fixed point impurity entanglement entropy is non-zero at the MG-point and attains 
its maximum of ln(2) in the middle of the chain. 




r/R 



Figure 14. DMRG results for S lmp {J' K = 0, r, R) for R = 100, 200 for the spin-chain 
model with J 2 = J/2. The solid line represents the theoretical result Eq. (|6.35l) . 

Likewise, we can obtain Si mp (J' K = 0,r,R) for R even at the MG-point using the 
TS-ansatz with variationally determined ipn° l - The results are indistinguishable from 
the OTS-ansatz, Eq. (16.351) . with p(r) from Eq. (16.381) . Our results obtained at the 
MG-point are shown in Fig. [T4]for system sizes R = 100, 200. Almost no variation with 
R is observed and the numerical DMRG results display almost complete agreement 
with the TS- and OTS-ansatz. The TS- and OTS-ansatz yields results that are almost 
identical and the difference between the two is not visible in Fig. [TH The lack of finite- 
size effects is presumably related to the fact that Si mp (J' K = 0,r,R) does not contain 
a contribution from the single particle entanglement entropy (SPE). Si mp (J' K = 0,r/R) 
display a distinct cross-over between ln(2) and zero and clearly remains non-zero as 
R — > oo. 

7.2. Fixed Point Entanglement in the Dimerized Phase, J 2 < J% < J~/2 

We now turn to a discussion of the variation of the fixed point entanglement entropies 
Sim P {.J'K = ^y r /R) an d S imp (J' K = 0,r/R) with J 2 as we decrease J 2 towards the 
critical point J|. In Fig. [15] we show DMRG results for S imp (J' K = l,r/R) for 
J2 = J^j 0-3 J, 0.4 J, 0.5 J. In all cases with R = 401. The dashed line represents the 
OTS-result, Eq. fl6\28l) . with p(r) from Eq. flOSl) . valid at the MG-point. It is seen 
that S imp (J' K = l,r/R) quickly diminishes as J 2 is decreased towards J|. At J 2 we have 
argued in section [5] that Si mp (J' K = 1, r/R) approaches zero in the thermodynamic limit. 
In fact, from the structure of the OTS-ansatz it is clear that this approximation must 
break down at the critical point since almost any form for the soliton wavefunction, i/j so1 , 
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Figure 15. DMRG results for Si mp (J' K = 1, r, R) for R = 401 for the spin-chain model 
with J 2 = J|,0.3J,0.4J,0.5J. The dashed line represents the OTS result Eq. 
valid at the MG-point (J 2 = J/2). 



would imply a non-zero p and thus a non-zero result in Eq. (16.281) . We also note that 
at the critical point the soliton becomes massless. 



ln(2) 



OTS-ansatz 

S imp (J 2 =0.5J,J' K =0,r,R=100) 
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Figure 16. DMRG results for Si mp (J' K = 0, r, R) for R = 400 for the spin-chain model 
with J 2 = J|,0.3J,0.4J,0.5J. The dashed line represents the OTS result Eq. (jQ5)) . 

We now focus on Si mp (J' K = 0,r/R) for which we show DMRG results for 
SimpiJ'x = l> r /-R) with J 2 = J|, 0.3 J, 0.4J, 0.5 J in Fig. [16j The dashed line represents 
the OTS-result Eq. (16.351) . with p{r) from Eq. (16.381) . valid at the MG-point. Although 
some variation of Si mp (J' K = 0,r/R) with J2 is seen, this variation is rather small 
and this fixed point impurity entanglement entropy clearly remains non-zero at the 
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critical point, J|. There are some finite size effects but at J 2 these affects are rather 
small. Perhaps surprisingly, the OTS result Eq. (16.351) is not dramatically different 
from the DMRG results even at the critical point J|. We speculate that this is due 
to the fact that S imp (J' K = 0,r/R) does not contain a contribution from the single 
particle entanglement but can be described entirely in terms of the impurity valence 
bond picture. S imp (J' K = l,r/R), on the contrary, is completely described by the single 
particle entanglement. 



7.3. Fixed Point Entanglement in the Gapless Heisenberg Phase, J 2 < J 2 




Figure 17. DMRG results for S lmp {J 2 = 0.2 J, J' K = 1, r, R) for R = 101, 201, 301, 401 
for the spin-chain model with J2 = 0.2J. The solid line is DMRG results for 
Simp{J L 2 — J^J'k = 1, .R = 401) calculated at the critical point, J 2 = J|. 

Finally we turn to the gapless Heisenberg phase where J 2 < Jf. Our results 
for S imp (J2 = 0.2J,J' K = l,r/R) are shown in Fig. [T7| for a series of system sizes, 
R = 101, 201, 301 and 401 all calculated at J' K = 0.2 J. For comparison we have included 
results for S imp (J2 = J^J'k = l> r /-R) & t the critical point. It is seen that S imp (J2 = 
0.2 J, J' K = l,r/R) is rather small and negative for this value of J2 with relatively little 
variation with the system size R. In comparison Si mp (J2 = J^J'k = ^-i r /R) remains 
positive. A likely explanation for this is that the single particle entanglement, largely 
characterizing Si mp (J' K = l,r/R), is very sensitive to J2 and no longer is a sensible 
quantity for J 2 < due to the gapless nature of this phase. 

Results for S irnp (J 2 = 0.2J, J' K = 0,r,R) for R = 100,200,300,400 are shown in 
Fig. HI] for J 2 = 0.2 J. We also show results for S imp (J' K — 0, r, R — 400) calculated at 
the critical point, J 2 = J 2 - I n this case the results at J 2 = 0.2J are very similar to the 
results at the critical point J 2 = J| and it seems possible that the results coincide in 
the thermodynamic limit. We speculate that the fixed point entanglement found at the 
critical point is representative for all of the gapless Heisenberg phase J 2 < J|. 
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Figure 18. DMRG results for S lmp { J 2 = 0.2 J, J' K = 0, r, R) for i? = 100, 200, 300, 400 
for the spin-chain model with J 2 = 0.2 J. The dashed line represents the OTS result 
Eq. (|Q5]) . The solid line are DMRG results for S imp (J 2 = J%,J' K = 0, r, R = 400) 
calculated at the critical point, Ji = J%- 



8. The alternating part, Sa 

8.1. General discussion 

Entanglement entropy in spin chains also has an interesting staggered part as we pointed 
out in [H7J . For an open chain with no impurity, this was shown to decay away from the 
boundary with a power law 1/2, the same power- law exhibited by the dimerization. The 
operator of dimension 1/2, representing the dimerization, is related by a chiral SU(2) 



transformation to the staggered spin density, n, reviewed in Appendix A It has no 
counterpart in a non-interacting fermion system. (Or, more correctly, the counterpart, 
while it exists, must always come together with the operator cos0 c from the charge 



sector as reviewed in Appendix A).) This alternating entropy for the open chain with 
no impurity decays with a different exponent than what occurs for an open chain of free 
fermions [87] . If we now include a weak coupling of an impurity spin to the end of the 
chain, we expect that the associated change in the alternating part of the entanglement 
entropy will continue to be different that of the free fermion Kondo model. This is in 
striking contrast to the uniform part which we have argued to be the same (at long 
length scales) for free fermion and spin chain Kondo models. In the rest of this section 
we discuss the behavior of the dimerization and examine the behavior of the alternating 
part of the entanglement entropy for both critical and dimerized spin chains. 

8.2. The Alternating Part of the Energy Density, Ea 



We start by deriving a field theory epxression for the alternating part of the energy 
density, Ea, that we shall find sheds some light on the alternating part of the 
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entanglement entropy. The energy density for XXZ antiferromagnetic spin chains: 

(h r ) = ((S+S; +1 + S~S+ +1 )/2 + AS z r S z r+1 ) (8.1) 

is uniform in periodic chains. On the other hand, an open end breaks translational 
invariance and there will be a slowly decaying alternating term or " dimerization" in the 
energy density 

(h r ) = E u (r) + (-l) r E A (r), (8.2) 

where E A (r) becomes nonzero near the boundary and decays slowly away from 
it. We can calculate E A (r) by Abelian bosonization modified by open boundary 
conditions [106] . In the critical region |A| < 1, the low energy effective Hamiltonian is 
just a free massless relativistic boson. 

The staggered part of h r ~ (-l) r+1 (^ L - ip[ip R ) ~ (-l) r+1 sin( v / 4^0). 
Here we follow the notation of Ref. [1U7] , but define the Luttinger parameter as 
K = 7r/(2(7r — cos -1 A)) so that K = 1 for an XY spin chain and K = 1/2 for 
the Heisenberg model. In a system with finite R and open boundary conditions, 

1 

tv^f)] J 



E A (r,R) cc (sin(v / 4^0)) cx ^— — — . (8.3) 



This is our basic result for E A , from which it follows that E A (r, R) = f{r/R)/yR for 
the Heisenberg model but E A (r,R) = f{r/R)/R for an XY spin chain, corresponding 
to free fermions. 

At the Heisenberg point, A = 1, Eq. (18.31) will have some logarithmical corrections 
due to the presence of a marginally irrelevant coupling constant, g, in the low energy 
Hamiltonian, Eq. flA.161) . that we now try to take into account. This interaction is 



reviewed in Appendix A Ignoring boundaries, the staggered energy density E A ~ 
sin(\/27r0) has the anomalous dimension 

y(g) = 1/2 - 3(7/4. (8.4) 

With a boundary, the renormalization group equation for E A is the naive one, involving 
the anomalous dimension, 7(g) in the usual way. 

[d/d(ki r) + /3(g)d/dg + y(g)]E A (r, r/R) = 0. (8.5) 

(Here the partial derivative with respect to r is taken with r/R held fixed.) 

This result may require some justification since, in general, the presence of a 
boundary can strongly affect the scaling behavior. It is crucial here that we are 
considering E A far from the boundary compared to the ultraviolet cut-off (i.e. r> 1). 
The boundary condition dictates that we should regard the right moving factor in the 
operator e ly fi^<l } ( r ) ^ which occurs here as a left moving operator at the reflected point — r: 

e iV2w(f>(r) _ e iV2^(<j} L (r)+cj> R (r)) _^ e iy/2^(<f> L (r)-<t> L (-r)) /g 

To calculate the anomalous dimension of this bi-local operator we can consider 
the operator product expansion (OPE) with the marginal interaction (reviewed in 



Appendix A). This marginal interaction, cos V87T0 also becomes bi-local in the presence 
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of the boundary. The OPE of these two bilocal operators, e iv/ ^^ i ^ r - ) ~* L< '~ r ^ and 
e -zV8n(<t>L(r)-<t>L(-r)) j g produce Q f the OPE's of the operators at ±r separately. 
Fortunately, this gives exactly the same result as without the boundary. Therefore 
we expect the naive RG equation to apply. 
The general solution of Eq. (18.51) is 

E A (r,r/R)=F[g(r),r/R}ex V {- / d(ln r'h^O]}, (8.7) 

Jro 

where F is an arbitrary function of g(r) and r/R. 
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Figure 19. Entanglement entropy S(r) (black squares) and alternating part of 
the energy density E A (r) (open circles) computed with DMRG at Jf for R = 200 
sites. The dashed line is Eq. (|8.3p . The uniform part of S(r), obtained by taking 
S(r) + (— l) r aEA{r) with a = 1.001699, is represented by open squares. The best fit, 
shown by a red curve, is indicated on the plot. 

In the weak coupling limit, we may evaluate the scaling function F(g(r),r / R) at 
g{r) = and with the second order beta function, (3 = —g 2 , we obtain: E A (r,r/R) = 
F{r/R)/ [t/t (In \r |) 3 ^ 4 ]. One can push this a bit further following a similar calculation 
in Ref . [98j 1108] . Provided with the beta function up to third order 

(3(g) = g 2 - (l/% 3 , (8.8) 

the effective coupling solved from Eq. ( 18.81) is 

{ln(r/r ) + \ m[m(r/r )]}, (8.9) 



g(r) 9o 2 
and expanding F(g(r),r / R) in powers of g(r), we can improve the solution as 

F(r/R) 



E A (r,r/R) 



y/r[\n(r/ai) + |lnln(r/ai)] 

(8.10) 



2 

« 2 



[ln(r/ai)p 
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where a term proportional to l/[ln(r/ai)] which could have occurred inside the curly 
brackets can always be adsorbed by redefining a\. Note that since r/R is being held 
fixed here, as we take r — > oo, it follows that we can always replace r by R inside the 
logarithms in Eq. (18.101) by rescaling a\ by R/r. 

As reviewed in Appendix A at the critical value of J2, the marginal coupling 
constant vanishes and all logarithmic corrections vanish. Hence, we can check Eq. ( 18.31) 
directly using our DMRG results if we work at the critical point J|. This is shown in 
Fig. [191 where results for Ea{t, R = 200) for a uniform (J' K = 1) system with R = 200 
are plotted along with Eq. (18.31) . The agreement is quite good. In Ref. [87] it was argued 
that the alternating part of the entanglement entropy is given by Sa = —oEa- With 
a = 1.001699 we show S(r, R) + (— l) r aEA{r, R) in Fig. [19] and excellent agreement with 
this quantity and the uniform part of S(r, R) is observed. Fitting this uniform part to 
Eq. (|2~3l) and using the fact that ^U7\ In g = -(1/4) In 2 we estimate si ~ 0.746. 



. 3. Scaling of at J, 
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Figure 20. (a) DMRG results for the total entanglement entropy, S(J' K ,r, R) for 
a 102 site spin chain at J|, with a J' K = 0.41 Kondo impurity (■) along with 
S(J' K = 1, 1 — 1, R— 1) (□). For both cases is the extracted alternating part shown along 
with the resulting Sf mp {J' K = 0.41, r, R = 102) for R even, (b) DMRG results for the 
total entanglement entropy, S(J' K , r, R) for a 101 site spin chain at Jf , with a J' K = 0.41 
Kondo impurity (■) along with S(J' K = 1, r — 1,R — 1) (□). For both cases is the 
extracted alternating part shown along with the resulting <S^ p ( J' K = 0.41, r,R = 101) 
for R odd. For comparison we also show S^AJ^ = 0.41, r 7 R~ 102) for R even from 
panel (a) (dashed line). 

Our fundamental definition of Si mp , Eq. ( 12.41) . focused only on the uniform part 
of the entanglement entropy. It is also possible to define the alternating part of the 
impurity entanglement entropy following Eq. ( 12.4ft : 

Sf mp (J' K , r, R) = S A {J' K , r, R) - S A {1, r-l,R-l). (8.11) 
As before, we have subtracted Sa when the impurity is absent, in which case both r 
and R are reduced by one and the coupling at the end of this reduced chain, linking 
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site 2 to 3 and 4, has unit strength. Applying this definition to numerical data involves 
some subtleties. First of all S A is only defined up to an overall sign. Secondly, when 
calculating Sa(1,t — 1,R — 1) we define this as — 5^(1, r', R') with R! — R — 1 since 
the shift from r to r — 1 implies a sign change in the alternating part. For convenience 
we have therefore always exploited this degree of freedom to use a sign convention that 
makes the resulting S^L, positive in all cases. In Fig. [20] we show data for the total 
entanglement entropy along with the extracted alternating parts and the resulting S^L, 
for both R = 102 even and R = 101 odd. The initial data are the same as shown in 
Fig. [3] with £ K = 25.65. As was the case for S imp we do not observe any special features 
in S^ p (r) for fixed R, J' K associated with the length scale and in all cases decays 
monotonically with r. 




r/4 K r/i; K 

Figure 21. y/rS^ np for fixed r/R and a range of couplings J' K at Jf. (a) for R even 
and (b) for R odd. The values of £k used are obtained from the scaling of Si mp and 
are listed in table [TJ DMRG results with m = 256 states. 

We now turn to a discussion of a possible scaling form for S^p. In Ref. [87] 
it was shown that for J' K = 1 the alternating part of the entanglement entropy, 
S A , is proportional to the alternating part in the energy, Ea and it was shown that 
Ea(t) = f{r/R)/yf{r) for some scaling function /. See also the detailed derivation 
in subsection 18.21 A first guess for a scaling for would then simply follow a 
generalization of the above formula to the case J' K ^ 1. Naively, this would imply 
that v^O^imp should be a scaling function, f(r/R,r/£ K ). Our results for v^O^imp 
for fixed r/R are shown in Fig. [21] for a range of J' K and R. The values for £k used 
to attempt the scaling are the ones previously determined from the scaling of Si mp for 
fixed r/R at J|, listed in tabled] Clearly the results for \frS A m p follow the expected 
scaling form. We expect that the scaling would have been better had we allowed the £k 
to vary instead of using the data from table [T] 
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Figure 22. DMRG results for the alternating part of the fixed point entanglement, 
S^ lp { J 2 , J' K = 1, r/R), shown for a range of values of J 2 = Jf , 0.3 J, 0.4J, 0.5 J, in each 
case for 2 different values of R = 101,401. The dashed line indicates the result for the 
OTS-ansatz, Eq. i[8T2|) . 



8.4- Alternating Part of the Fixed Point Entanglement Entropy 

Following our definition of the alternating part of the impurity entanglement entropy, 
S^L,, given in Eq. (I8.1ip . we can also analyze the alternating part of the fixed point 
entropy by studying Sf mp {J' K = 0) and Sf mp (J' K = 1). If we first focus on Sf mp (J' K = 0) 
it is clearly zero when R is odd as was the case for Si mp (J' K = 0). For R even, we 
see from Eq. flfTMD and Eq. fl6\27D that at the MG-point (J 2 = J/2) the OTS-ansatz 
predicts that also in this case is 5*^ p (J^ = 0) zero. Numerical DMRG results confirms 
this and shows that also for J 2 < J/2 is (J' K = 0) negligible. We have verified that 
this also holds at the critical point J2 = J| as well as in the gapless Heisenberg phase 
at J 2 = 0.2J. 

The interesting fixed point entropy is then S^ ip (J' K = 1). As for Si mp {J' K = 1), 
we see that the difference between R even and R odd amounts to a change of sign of 
the resulting fixed point entanglement entropy. We therefore only consider R odd. At 
the MG-point we expect the OTS-ansatz to yield very precise results. If we combine 
Eq. (I6.27P with the fact that for R even S A (J' K = 1) = ln(2)/2 for any r, we find the 
OTS result: 

SLM'k = 1) = (1 - P) ln(2), R odd. (8.12) 

Where p is given in Eq. ( 16.381) . As we mentioned in the previous section there is a 
slight subtlety here since, in order to calculate S^^J'K = 1) for R odd, we need 
5^(1^ — 1, R — 1) for R — 1 even. Due to the shift, r — * r — 1, we have defined 
this as — ln(2)/2 as opposed to the ln(2)/2 we just quoted. This sign convention 
renders the resulting fixed point entanglement entropy positive for all r/R. Choosing 
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the other possible sign convention would have resulted in a trivial shift in the fixed point 
entanglement entropy of hi(2), yielding — phx(2). 

Our numerical DMRG results for S^^J^ = 1) for R odd as a function of r/R are 
shown in Fig.[22]for a range of second nearest neighbor couplings J 2 = J%, 0.3 J, 0.4 J, 0.5 J 
throughout the dimerized phase J| < J 2 < J/2. In each case we show results for two 
system sizes R = 101 and 401. At the MG-point we see that finite size effects are minimal 
and the results for R = 101 and 401 agree very well. In both cases the agreement with 
the result from the OTS-ansatz, Eq. (I8.12p . shown as the dashed line, is almost perfect. 
As J 2 is decreased from the value of J/2 pronounced finite-size effects develop and the 
fixed point entanglement entropy clearly tends toward zero with increasing R. This 
is also true at the critical point, J 2 = J 2- F° r clarity we have not shown results for 
J2 = 0.2 J in the gapless phase where (J' K = l,r/R) for both R = 101 and 401 is 
smaller than at the critical point. 

9. Conclusions 

The presence of impurities will clearly affect the entanglement and we have here defined 
the impurity contribution to the entanglement entropy, Si mp . Using the equivalence of 
the electronic Kondo model and the J\ — J2 spin chain model in the regime J2 < J% 
we have shown numerical evidence that S imp follows a scaling form S imp (r /£ K ,r / R) 
demonstrating the presence of the length scale £k associated with screening of the 
impurity. We have provided rather strong arguments in favor of this scaling and 
analytical results based on a Fermi liquid approach valid for r ^> The D-dimensional 
Kondo model has been shown to yield the same entanglement as the 1-D Kondo model 
which in turn is equivalent to the entanglement in the J\ — J2 spin chain. The single-site 
entanglement obtained by only considering the entanglement of the impurity spin has 
been shown to display weak scaling violations. For J 2 > J% in the dimerized phase we 
have shown that almost exact calculations are possible using the variational TS-ansatz. 
Simplifying the TS-ansatz we have shown that the contributions arising from impurity 
valence bonds (IVB) and single particle entanglement (SPE) naturally arise at the MG- 
point. Finally, we have argued that at finite temperatures the entanglement entropy for 
a sub-system A of size r, S(T), will approach the thermal entropy, S t h(T), for T ^> v/r. 
In light of the generality of our model we expect our results to be quite widely applicable. 
The study of quantum impurity entanglement as it occurs in more complex models with 
many body ground-states displaying non-trivial order would clearly be of considerable 
interest . 
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Appendix A. Review of field theory results on the Kondo effect and spin 
chains 

The usual Kondo model describes electrons moving in 3 dimensions with a short range 
exchange interaction with a single S=l/2 impurity spin. No other electron-electron 
interactions are taken into account. Although similar conclusions can be arrived at 
more generally, the simplest case is where a free electron quadratic dispersion relation 
is assumed, and the Kondo interaction is taken to be a (^-function, implying spherical 
symmetry. The Hamiltonian is given in Eq. (11.31) . 

Kondo model calculations are traditionally carried out by reducing the 3D fermionic 
model of Eq. ( 11. 3D to an equivalent ID fermionic model. This can be done by expanding 

^(f) = (1/r) XXmWUr), (A.l) 

where the Y^ m (f) are the usual spherical harmonics, depending only on the direction off* 
and the ipi im (r) are ID quantum fields. The spherically symmetric H of Eq. (11.31) reduces 
to a sum of commuting terms: H = Ylim-Hi,m- Due to the 5-function interaction, all 
harmonics are non-interacting except for the s-wave. The higher harmonics contain only 
the centrifugal potential energy VJ(r) = 1(1 + l)/2r 2 . A low energy form of the / = 
Hamiltonian is obtained by linearizing the dispersion relation around the Fermi surface, 
e(k) = k 2 /2m — e_p, yielding a ID Dirac fermion theory on the half-line, r > 0, coupled 
to the impurity spin at r = 0. We introduce left and right movers by: 

^ =0 (r) « e lkFr ij R (r) + e^^r), (A.2) 

where if)L/n{ r ) vary slowly and /cj? is the Fermi momentum. 

t-R 



H XD « (iv/2ir) I dr[^ L {d/dr)^ L - ^ R (d/dr)^ 
Jo 

+ v\ K ii(0)(a/2)i> L (0)-S. (A.3) 

Note that we have adopted an unconventional normalization for these ID fermion fields: 

{^ L/R (r)^ L/R (r')} = 2n5(r-r'). (A.4) 

Here Xk oc Jk and the left and right movers obey a boundary condition: "0l(O) = ^(0) 
[771 1781 [79] . Up to some trivial geometric factors of r, all the physics of Eq. ( 11.31) can 
then be obtained from this ID fermion model. If the 3D system is originally defined 
inside a sphere of radius R, then the ID system exists on a line of length R. 

It is now very convenient to bosonize this ID model. This allows for the introduction 
of separate bosonic fields representing the spin and charge degrees of freedom of the 
conduction electrons. The non-interacting free electron Hamiltonian can be written as 
a sum of decoupled spin and charge terms: 

H = H s0 + H c0 . (A.5) 
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The spin part can be written in terms of the spin current operators: 

Jl/R = ^ L/R ^L/R, (A.6) 



as: 

-R 



H s0 = (v/2ir) [ dr(l/3)[J L ■ J L + J R ■ J R \. (A.7) 
J o 

The free boundary condition on the fermion fields at r = implies a boundary condition 
in the continuum limit: V'l(O) = — V'-r(O) which in turn implies Jl(0) = Jr(0). Since 
the Kondo interaction only involves Jl it can be written entirely in terms of the spin 
boson field. 

H K = v\ K J L (0) ■ S. (A.8) 

(Actually, there are irrelevant operators which couple spin and charge sectors together, 
but these can be ignored in the low energy effective theory.) 

It is actually possible to continue to the negative r axis, defining: 

M-r) = -Mr), (0<r<R). (A.9) 
This in turn implies that Jr(t) = Jl{— t). The free spin Hamiltonian then becomes: 

H s0 =(v/6n) [ drJ L -J L . (A. 10) 

J-R 

Imposing convenient boundary conditions on left and right fields at r = R, we obtain 
the periodic boundary conditions in the left-moving formulation: Jl(—R) = Jl(R). 

Now consider the J\-Ji spin chain model of Eq. (11.41) . The connection between 
spin chains and the Kondo model was pointed out in [1071 I109j . using field theory 
arguments and in |100j the Hamiltonian of Eq. (11.41) with J 2 = was, among others, 
solved exactly by Bethe ansatz noting connections with Kondo physics. The low energy 
theory can again be described by bosonization. One approach is to bosonize the weakly 
coupled Hubbard model at half-filling and then to extrapolate to strong coupling. We 
again obtain spin and charge bosons. The Hubbard interaction gives a gap to the charge 
boson field, which can thus be eliminated from the low energy theory. The spin operators 
at site j can be represented: 

Sj ~ [Jl + Jr] + (— Inconstant cos (j) c ■ n. (A. 11) 

Here C is the charge boson and n is the antiferromagnetic order parameter and can 
be written entirely in terms of the spin boson field. The Hubbard interaction makes 
< cos 0c >t^ so we may replace this factor by a constant: 

Sj w [J L + J R ] + (-Inconstant • n. (A. 12) 

The low energy Hamiltonian is simply the spin part of the non-interacting electron 
Hamiltonian, H s0 , up to irrelevant operators. In particular, the spin current operators 
have the same Green's functions as in the non-interacting model. The field n can also 
be expressed in terms of the non-interacting spin bosons. It is represented in terms of 
exponentials of bosons (in the usual abelian bosonization scheme). 
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Now suppose that the chain has a free end at r = 0. The boundary conditions on 
the fermion fields can be translated into boundary conditions on the spin boson field. 
These can be seen to imply: 

n(0) oc J L (0) = 4(0). (A.13) 

Now suppose that we weakly couple one additional spin to the spin chain, with 
a coupling constant J' k <C 1. We may write the low energy effective Hamiltonian 
using bosonization. If the extra spin couples at the end of the chain then the resulting 
interaction term is particularly simple: 

H K oc J L (0) • S. (A. 14) 

The coupling constant here is proportional to J' K but the constant of proportionality 
is non-trivial involving the proportionality constant in Eq. (1A.13I) . It can be extracted 
by studying the end-to-end spin correlation function in the chain without the weakly 
coupled spin [55], yielding: 

J' K « 1.3807A*. (A.15) 

Thus the same low energy effective Hamiltonian is obtained for the spin chain with a 
weakly coupled spin at the end as for the usual free fermion Kondo model. 

Two caveats should be made here. First of all, if the weakly coupled spin couples 
to a point far from the end of the spin chain then a much different low energy effective 
Hamiltonian occurs. Secondly, even when the weakly coupled spin is at the end, there 
are subtle differences from the usual Kondo model due to differences in the irrelevant 
operators. In the free fermion model all irrelevant operators have dimensions of at least 
3 and are strictly irrelevant. However, the spin chain model has a notorious marginally 
irrelevant operator in its low energy effective Hamiltonian. 

SH = -(gv/27r) J L ■ J R . (A. 16) 

This leads to logarithmic corrections to essentially all low energy properties which 
are difficult to calculate in detail and cause problems with fitting numerical data. 
Fortunately, there is a way of avoiding this difficulty. As we vary J 2 , this marginal 
coupling constant, g also varies. For small J 2 , g > 0, and is marginally irrelevant. A 
critical point occurs at J 2 = J c « .2412. For larger J 2 the system goes into a gapped 
spontaneously dimerized phase. Right at the critical point g = and all logarithmic 
corrections vanish. Much better agreement between numerical simulations and field 
theory predictions can be obtained at this point. All remaining irrelevant operators 
are strictly irrelevant, leading only to corrections which vanish with power-laws of the 
energy (or inverse length) scale, not logarithms. In a separate paper [55], we study the 
J 2 = Kondo spin chain model in more detail. Here we will focus on the J 2 = J c case. 
Note that this is actually much closer to the free fermion version of the Kondo model. 
We note that Eq. flA~T5l) was obtained at J 2 = J c . 

At low energies and long length scales, the effective Kondo coupling becomes 
large, and the effective Hamiltonian flows to the strong coupling fixed point. In the 
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electron version of the Kondo model we may think of this fixed point as one where 
the impurity spin is "screened", i.e. it forms a singlet with a conduction electron. 
The remaining electrons behave, at low energies and long length scales, as if they 
were non-interacting, except that they obey a modified boundary condition reflecting 
the fact that they cannot break up the singlet by occupying the same orbital as the 
screening electron. This modified boundary condition corresponds to a ir/2 phase shift. 
Correspondingly in the spin chain Kondo model, the impurity spin gets "adsorbed 
into the chain" and no longer behaves like a paramagnetic spin at low energies and 
long distances. The leading corrections to this low energy long distance picture are 
described by lowest order perturbation theory in the leading irrelevant operator at the 
strong coupling fixed point. This is an interaction between the remaining conduction 
electrons, near the screened impurity. (It doesn't involve the impurity itself since it 
is screened and doesn't appear in the low energy effective Hamiltonian.) This leading 
irrelevant operator is J L (0) • Jr. (0) [781 175]. From flXTTDft . we see that this is proportional 
to the spin part of the free electron energy density. It is the entire energy density 
in the low energy effective Hamiltonian for the spin chain. The energy density has 
dimensions of (energy) / (length) so the corresponding coupling constant in the effective 
Hamiltonian must have dimensions of length. On general scaling grounds we expect it 
to be proportional to The precise constant of proportionality simple corresponds to 
giving a precise definition of what we mean by We adopt the convention: 

H int = -«k)H s , l (0). (A.17) 

Here the subscripts s and L are a reminder that this is the spin only part of the energy 
density for left movers. For the purpose of doing first order perturbation theory in 
H int f° r quantities like the susceptibility, specific heat or ground state energy, which are 
translationally invariant in th order, we may replace [77J Hi nt by: 

H^t -Hk/{2L)\ I H s , L {x). (A.18) 

J-R 

This is equivalent to a length dependent reduction of the velocity: 

v^v[l-n^ K /(2R)]. (A.19) 

This then implies that the susceptibility, which is Rj {2nv) in the absence of the Kondo 
impurity becomes: 

X - R/{(2irv)[l - ^ K /(2R)]} « R/(2ttv) + £ K /(4v) 

= R/(2nv) + 1/{AT K ). (A.20) 

Thus the zero temperature impurity susceptibility is 1/(47V). It is this form of the 
impurity susceptibility, simply related to the high temperature, free spin behavior, 
1/(4T), which motivates the definition of (and hence T K = Vf/£k) implied by 
( 1A.17I) . We note that this interaction Hi nt is present even in the absence of an impurity, 
for free fermions but then the coupling constant is of order a lattice constant. Similarly 
it is also present for the spin chain with no impurity (i.e. J' K — 1) with a coupling 
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constant of order a lattice constant. The effect of a weak Kondo coupling is to make 
this coupling constant large. We emphasize that this precise choice of definition of Tk 
has no physical consequences. The power of Fermi liquid theory is to predict not only 
the form of low energy quantities but also ratios of various low energy quantities such as 
impurity susceptibility, impurity specific heat, resistivity, etc., corresponding to various 
generalized Wilson ratios. In the limit £g- <C r, we can also calculate Si mp using the 
FLT interaction of Eq. IA. 171 in the lowest order perturbation theory. 

Appendix B. 3D to ID reduction and entanglement entropy 

In this appendix we prove that the impurity entanglement entropy is the same for the D 
and one-dimensional Kondo model and consequently also for the spin chain model. We 
also present a new derivation of the free fermion entanglement entropy in D-dimensions. 

With a spherically symmetric dispersion relation and Kondo interaction, as in Eq. 
(11. 3p . the Hamiltonian is a sum of terms acting on subspaces of different I, m quantum 
numbers for single electrons. That is, the kinetic energy term, for each electron separates 
into such a sum and so does the Kondo interaction. If the Kondo interaction is a 8- 
function then only the s-wave term is non-zero. If it is longer range there are small 
larger / terms. However, according to the usually RG picture of the multi-channel 
Kondo model the largest term grows under renormalization and the others shrink, so 
that in the low energy effective Hamiltonian we can ignore the higher harmonics. 

This form of H implies that the ground state wave-function is a product of factors, 
one from each channel. This follows because the single-particle wave-functions all have 
definite I, m quantum numbers. We may define electron creation operators c\ mr which 
create an electron at distance r from the origin in the I, m channel, and write the ground 
state as a product over I and m of factors for each /, m. Ignoring the Kondo interaction 
(as we can basically do for I > 0) each factor just corresponds to the product of creation 
operators corresponding to a filled Fermi sea. (Each channel is filled up to the Fermi 
energy.) Including the Kondo interaction this factorization is still valid but the wave- 
function for channels with a Kondo interaction are non-trivial. So, of course, the pure 
density matrix \tp >< ip\ is also a product of factors for each /, m. 

Now consider tracing over the region outside a sphere of radius r. It is important 
here that we choose a sphere, not some other shape, so as to preserve the rotational 
symmetry. This should leave a reduced density matrix which is also a product over / and 
m. This may be especially clear if we introduce an ultra-violet cut-off by only allowing 
the particles to sit on the surfaces of spheres at various distances from the origin. The 
reduced density matrix is obtained by tracing over all the spherical surfaces further from 
the origin than r. So we can write: 



Here each reduced density matrix, p^ m acts on a different sector of the Hilbert space. 
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So we can write: 

lnp(r) =^lnp l>Tn (r). (B.2) 

Z,m 

Actually, this notation implies that each term is a product of a non-trivial operator in 
one channel and the identity operator in all the others. When we calculate trp In p, we 
get a sum of terms like trp In p/ >m . Because we can write a complete set of states as a set 
of products in each channel, the trace of an operator written in product form reduces 
to a product of traces. So: 

Trplnp /jm =[ Trp^ m /] Tr(p ljm lap l>m ). (B.3) 

Since trp^ m / = 1, we get: 

Tr[plnp] = y^Tr[p; >m lnp; im ]. (B.4) 

Now if we compare the zero Kondo coupling to finite Kondo coupling case, only 
the s-wave channel changes. (More correctly, for a finite range Kondo interaction there 
is also some change in the higher / channels. However, this is presumably just a short 
distance effect and drops off much more quickly than the s-wave part. i.e. like a/r 
where r is a cut off scale, instead of £,k/t for the s-wave.) So the impurity entanglement 
entropy for the 3D system is the same as for ID. Exploiting the equivalence of the ID 
Kondo model and the spin chain model the entanglement entropy is the same also for 
this model. 

We remark that the well-known r 2 In r form of the entanglement entropy [1101 IX 111 
I112j for free fermions (at R — > oo) can apparently be recovered from the decomposition 
into angular momentum channels. Note that the Hamiltonian for all / > contains a 
centrifugal potential: 

H ^ = -2d^ + ^P^- (R5) 
This centrifugal potential will not be too important when the size, r, of region A is 
sufficiently large and / is not too large. Thus we expect to recover essentially the usual 
ID result S = (1/6) In r for these values of / and m. However, for large enough /, we 
expect the entanglement entropy to be reduced, and eventually to vanish at large / for 
any fixed r. This can be seen from the fact that the density of electrons in region A 
with angular momentum I vanishes at large I. If the region A is essentially empty there 
can't be any entanglement. We may estimate the order of magnitude of I at which the 
sphere of size r becomes empty by simply comparing the centrifugal potential to the 
Fermi energy, 

,-, lmax\lmax 1) /t-, 

Ef ~ 4^ • ( R6 ) 

This gives 

l max w 1r^E~ F . (B.7) 
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Thus we estimate the free fermion entanglement entropy in 3D as: 

Imax 771 2 

S(r) « (1/6) lnr V(2Z + 1) « lnr. (B.8) 

3 

£=0 

(The prefactor of -Ep oc n -2 / 3 where n is the electron density is only expected to be 
correct in order of magnitude. ) This argument generalizes to any dimension D, giving: 

Soar - 1 lnr. (B.9) 
Appendix C. Separating Staggered and Uniform Parts 

In this appendix we derive the 7-point formula used for extracting the uniform and 
alternating parts of the entanglement entropy from the numerical data. 

We focus on functions, /, defined on a set of discrete lattice sites i, with a uniform 
part, u, in addition to a staggered part, s: 

/(<) = «(<) + (-i)«*(0. (c.i) 

We assume that both s and u are slowly varying. Often s and u are extracted by using 
the simplest possible 2-point approximation 

/(0 + /(< + !) 



WU 



2 

s(z) = (-1) , (C.2) 

by effectively assuming that u(i) ~ u(i + 1) and s(i) ~ s(i + 1). This somewhat 
crude approximation has many drawbacks and is insufficient for the present study. We 
therefore focus on higher order n-point approximations by approximating u and s locally 
by polynomials. If one is interested in developing an n-point formula that is symmetric 
around the point of interest it is not possible to use the same degree of polynomial for 
both u and s. For the present study the most frequently occurring case is u varying more 
rapidly than s. We have then found it sufficient to develop a 7-point approximation by 
assuming that u(i) ~ ai 3 + bi 2 + ci + d and s(i) = ei 2 + fi + g, arriving at the equations: 

f(i - 3) ~ - 27a + % - 3c + d - (9e - 3/ + g) 
f{i - 2) ~ - 8a + Ab - 2c + d + (4e -2f + g) 
f(i - 1) ~ - a + b - c + d - (e - / + g) 
f{i) ~d + g 

f(i + l)~a + b + c + d-(e + f + g) 

f(i + 2) ~ 8a + Ab + 2c + d + (4e + 2f + g) 

f(i + 3) ~ 27a + 9b + 3c + d - (9e + 3/ + g) 

(C.3) 
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Solving these equations for d, g we immediately get: 

u(i) = - —f(i - 3) - —f(i - 2) + —f(i - 1) 

1137 i i 

+ i /W + ^/( i+ D + ^ + 2) --/(, + 3) 

1137 1 1 

+ LfU) _ —f(i + i) /(i + 2) + — f(i + 3)1 

T v ; 496 ,/ v ; 248 ,/ v ; 3r v n 

= (-l)«(/(t) - «(<)) (C4) 

Appendix D. Detailed Proof of limited SU(2) invariance of S 1 

In this appendix we show that the von Neumann entropy for a system in doublet state 
(R odd) is SU{2) invariant. For R even the ground-state is a non-degenerate singlet an 
S is manifestly SU(2) invariant. 

Consider the case of R odd, where the ground states form a spin doublet. The most 
general linear combination of the doublet of ground states can be written as a unitary 
transformation of the spin-up ground state: 

\rl> >= U\ |> • (D.l) 

Here we can choose 

U = exp[i(a x S% + a y S y T )} (D.2) 

so that 

U\ t> = [cos|a/2|Z+2S-£ T sin|a/2|]| t> 

= cos|a/2|| t> + (ax ~, my) sin|a/2|| |>, (D.3) 

\a\ 

with I the identity matrix and 

R 

and we have used the fact that the states | |>, | |> transform like an S=l/2 doublet 
under SU(2) to treat the Taylor expansion of the exponential. Now we observe that: 

U = U A ® U B , (D.5) 

where 

r R 

Ua = exp[m ■ Sj], Ub = exp[za • Sj]. (D.6) 

j=l j=r+l 

Thus: 

\^>=U A ®U B \]> ■ (D.7) 
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The pure density matrix can then be written: 

p = U A ®U B \]><]\U B ®U ] A . (D.8) 
It follows that the reduced density matrix is: 

p u A = U APA U\ (D.9) 



where p A is the reduced density matrix obtained from the pure state | |> and p A is the 
reduced density matrix for the unitarily transformed state U\ ]>. To prove Eq. (ID. 91) . 
consider decomposing the pure state | f> into a sum of a product of complete bases of 
states in regions A and B: 

| ]>=^C ab \a> ®\b> . (D.10) 

Then the pure state density matrix is: 

Pu=J2 CabC:, b ,U A \a > U B \b >< b'\U B < a'\U\ (D.ll) 

aba'b' 

To obtain the reduced density matrix we perform the partial trace over region B: 

p u A = C ab C* a , b ,U A \a><n\U B \b><b'\U B \n><a'\U\ (D.12) 

aba'b'n 

Here Ylm ls a sum over a complete set of states in region B. If we now use the following 
identity: 

< n\U B \b >< b'\U B \n >= ixU B \b >< b'\U ] B = tr\b >< b'\, (D.13) 

n 

we see that: 

Pa = U a J2 C abC* a , b ,\a >< a'\tv(\b >< b'\)U ] A = U A p A U ] (D.14) 

aba'b' 

Eq. ( 1D.9f) implies that p u A and p A have the same eigenvalues, and hence the same 
entanglement entropy. 

Appendix E. Connection between S and the thermodynamic entropy 

In this appendix we give arguments and provide numerical data showing that 
quite generally, for sufficiently high T, the entanglement entropy will approach the 
thermodynamic entropy. 

As usual, let us consider a subsystem, A of linear extent r with the entire system, 
A + B, of size R. Let us also consider the subsystem A decoupled from B. Following 
Gibbs, it is well known that the thermal entropy for subsystem A (decoupled from B) 
can be expressed as Sth = — Tr[pln p] where p = e~ HA ^ T /Z is the thermal density matrix 
and H A is the hamiltonian describing A. This mixed state is expected to occur if we 
start with our system A weakly coupled to an infinite reservoir, and then trace over the 
reservoir. Before performing the trace, we could regard the system as being in a pure 
state of system plus reservoir. On the other hand, the finite T entanglement entropy, 




Figure El. The finite temperature entanglement entropy per unit length of the 
subsystem for an XX spin chain shown along with the thermal entropy for system 
sizes of r = 2, 20, 50, 98. For the entanglement entropy the total system size was 
R = 100. All entropies are plotted per unit length. 



S = — tr[pA^- Pa\ is defined by beginning with the system A + B in the mixed state 
corresponding to the thermal density matrix e~( HA+HB >' T / Z and then tracing over region 
B to obtain the reduced density matrix p^(T). Again, we could arrive at this thermal 
density matrix by beginning with sytem A + B weakly coupled to an infinite reservoir 
and then tracing over the reservoir. While S t h{T) is, by construction, independent of the 
coupling to B, the entanglement entropy S(T) can clearly depend on it. We shall argue 
that for T ^> v/r the coupling to B can also be neglected when calculating S{T) and 
the entanglement entropy S will then approach S t h- Such a connection was previously 
noted in [26] and [ITS] . 

Essentially, we want to argue that we may regard region B approximately as a kind 
of additional reservoir for region A when T is sufficiently large. Normally a reservoir 
should be sufficiently weakly coupled to the system so as not to disturb its energy 
eigenvalues or eigenstates. This is not true in the case of subsystem B which is strongly 
coupled to subsystem A. However, this strong coupling only exists at the boundary 
(at the point r). For most eigenfunctions we can neglect the perturbation due to this 
coupling. The only important exception to this statement are the low energy states, with 
wave-lengths of order the system size and energies of 0(v/r). We can expect these to be 
strongly perturbed by the coupling to subsystem B and consequently expect S(T) to be 
strongly modified from the thermal entropy, Sth{T) at low temperatures. However, when 
T ^ v/r, we expect these states to make a negligible contribution to S(T) since the 
density of states is much larger at higher energies. Therefore we expect S(T) to approach 
Sth(T) when T ^> v/r. This argument suggests that this should occur regardless of R. 
We also expect this correspondance to hold regardless of boundary conditions. The 




Figure E2. The finite temperature entanglement entropy for subsystems r = 2,98 
for a R = 100 site XX spin chain. At T = the two entanglement entropies are 
identical and the dashed line indicates the difference between the two entropies, quickly 
approaching as T — > 0. 



thermal entropy for region A becomes independent of whether it is calculated with pbc 
or obc when T ^> v/r. Similarly the entanglement entropy becomes independent of 
whether region A is part of a system, A + B which obeys pbc or obc at T ^> v/r. 

In order to provide numerical evidence that S(T) indeed does approach S t h, the 
thermal entropy, we have performed exact calculations for both quantities on an XX 
spin chain of length R = 100. The finite T entanglement entropy is calculated for a 
subsystem, A, of size r within a total system A + B of size R obeying periodic boundary 
conditions. We have also calculated the thermal entropy, S t h(T), for a system of size 
r also with periodic boundary conditions. We then compare S(T)/r to S t h(T)/r. Our 
results are shown in Fig. IE1I where both entropies are plotted per unit length of the 
subsystem. Four different sub-system sizes of r = 2, 20, 50, 98 are considered and in all 
cases do we observe excellent agreement with the thermal entropy at sufficiently high 
temperature. In Fig. IE2I is shown the finite temperature entanglement entropy for a 
R = 100 site XX spin chain for two different subsystem sizes of r = 2, 98. At T = the 
two entanglement entropies are identical and the difference between the two is shown 
as the solid line, quickly approaching as T — > 0. 

Appendix F. Solitons in the Majumdar-Ghosh Model 

In this appendix we focus on various properties of the MG model that also can be 
calculated with very high precision using the TS-ansatz as described in section 16.11 We 
begin by a determination of 
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Appendix F.l. Determination of ip™ 1 

It is important to note that even though the thin soliton states \n) are linearly 
independent they are not orthogonal. The matrix of overlaps, B, is given by: 

B n<m =(n\m) = 2-\ n - m \. (F.l) 

If we now consider the action of the Hamiltonian Eq. (11.41) within the subspace spanned 
by the TS-states and define H' = H + 3JR/8 I (I is the identity matrix) it can be 
shown that [T03lfT04] : 



H'\0) «^[2|0>-|1>], n = 



H'\n) « - 
1 1 4 



i \ l \ i \ 

-\n) — \n — 1} — \n + 1} 



n + 0, N d 



H'\N d )^j[2\N d )-\N d -l)], n = N d . (F.2) 

Note that H' applied to a TS-state will in general also generate a contribution that is 
not in the TS-subspace. Here we have dropped the part of such "fat soliton" states 
that is orthogonal to the TS-states since we have restricted H' to the TS-subspace. We 
rewrite, Eq. ( 1F.2I) in the following way: 

H'\n) = jKJm), (F.3) 

where h' nm is a symmetric matrix with h' 0Q = h' N N = 2 and h' nn = 5/2 for n ^ 0, N d and 
h'nn+i = h'n+i n = — 1 for n = 1, . . . , — 1. Using the full hamiltonian H in Eq (11.41) 
we can then write: 

(*ls\H\*ls) _ ^=0 r^Am^f 3 

where, perhaps surprisingly, it turns out that the A m ^ n is simply given in matrix 
notation by A = h'BJ/4 = 3J//8. We can now determine the in a variational 
manner by minimizing the energy, Eq. ( IF. 41) . This is straightforward to do using the 
method of Lagrange multipliers and the constraint J2 m n '■Bm.nC 1 = 1, coming from 
normalizing the wave-function. Due to the simple form of the matrix A and the fact that 
the matrix B has only positive definite eigenvalues, it is then seen that the optimal values 
for the ^ s ° l is given by the eigenvector of B corresponding to the largest eigenvalue, X ma x- 
Furthermore, since {H') coincides with our definition of the soliton mass, we see that 
for open boundary conditions (obc) A^' obc is given by: 

Al s r = ^. (F.5) 

For R = 201,401 we find A^' obc /J = 0.12523,0.12506, respectively, in excellent 
agreement with the result for periodic boundary conditions. 

We note that the same ip^ oL can be found by assuming that (n\m) = S n>m and then 
solving the tight binding model resulting from Eq. (IF. 21) . In matrix form this tight 
binding model will simply be given by the matrix h'. However, since h'B is proportional 
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Figure Fl. Variationally determined ipn° l (open circles) for R = 201 shown with the 
the free particle form given by Eq. (|6.5p . 



to / an eigenstate of h! is also an eigenstate of B^ 1 {B^ 1 exists since B has positive 
definite eigenvalues) and then also of B. Hence, finding the eigenvector with lowest 
eigenvalue of h! will also find the eigenvector with the largest eigenvalue of B and hence 
the optimal if)™ 1 . This somewhat surprising observation, implies that at the MG-point 
the non-orthogonality of the thin soliton states cannot play an important role as we 
shall discuss in more detail later. 

As an illustration we show in Fig. IFll the variationally determined for R = 201 
along with the free particle form, Eq. (16. 5ft . Clearly the agreement is relatively good 
even for this relatively small system size. 

Appendix F.2. S* for R odd (J' K = 1) 

We now turn to a discussion of (S*) which also can be calculated employing the TS- 
ansatz, showing several surprising features. 

The TS-ansatz, Eq. (16.21) . is constructed using the single soliton states which do 
not include states where the soliton is on an even site when R is odd. It would then 
seem natural to assume that if we calculate (S*) using the TS-ansatz it would then be 
zero for r even. However, due to the non-orthogonality of the thin soliton states this 
turns out not to be the case. In fact we find that: 




(F.6) 



with, for r odd: 




otherwise 



(F.7) 
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Figure F2. (S?) at the MG point ( J 2 = J/2) calculated using the TS-ansatz, Eq. ([Oj) 
and variationally determined ifjn° l f° r ^ = 201 (solid circles). DMRG results for (S*) 
at the MG point keeping m = 256 states (open circles). 

and for r even 

{_l 2 -\n-m\ n< ^ m >r_ 
_l 2 -\n-m\ m< r_^ n yr_ (Rg) 
otherwise 

In Fig. [F2] we show results for (S*) at the MG-point (J 2 = J/2) obtained using 
the TS-ansatz as well as numerically using DMRG methods. Fairly good agreement is 
observed. Notably, the TS-ansatz clearly yields a non-zero value for (S^) when r is even. 
However, longer valence bonds, not accounted for in the TS-ansatz, clearly contribute 
for this value of R. 

We also stress that, if we use the OTS-ansatz to describe (S*) we simply find 
{S*) = \fp^°!_iy 2 \ 2 f° r r f° r r even. The non-zero result for for even r obtained 
with the TS-ansatz (See Fig. IF2[) is then purely a result of the non-orthogonality of the 
TS-states. 



Appendix F.3. Dimerization for R-odd from the TS-ansatz 

Finally, at the MG-point, J 2 = J/2, it is also possible to obtain very precise results for 
the dimerization using the TS-ansatz. Most measures of the dimerization are based on 
the spin correlation function. We therefore consider calculations of (S r ■ S r+ i) using the 
TS-ansatz. With r = 2n + 1, n = . . . N d = (R- l)/2 and B n , m = 2^ n ~ m ^ we find that: 

X^ N d J, sol p I sol 

. = y^p^ , (f.9) 
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Figure F3. (S r ■ S r +i) at the MG point (J 2 = J/2) calculated using the TS-ansatz, 
Eq. (|6.2p and variationally determined ?/>„ for i? = 201 (solid circles). DMRG results 
for (5 r 5 r +i) at the MG point keeping m = 256 states (open circles). 



with, for r odd: 



_3o-|n-m| „ > r^l or m > 1=1 

otherwise 
and for r even 

F nm = { -I 2 "'""" 1 ' ™<§orm<§ _ 
otherwise 

In Fig. [F3] we show results for (S r ■ S r +i) at the MG-point (J 2 = J/2) obtained 
using the TS-ansatz as well as numerically using the DMRG methods. Almost perfect 

agreement is observed. 
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